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Abstract 

A  given  report  is  based  on  the  results  of  studying  a  stationary  plasma-low  power-thruster 
(SPT)  and  an  a  ablative  pulsed  plasma  one  (APPT)  made  in  accordance  with  the  calendary  plan  of 
the  studies  on  the  ISTC-Project  No.23 1 7p  for  the  period  from  April,  1 ,  2002,  to  March,  3 1 , 2003 . 

In  the  first  part  of  the  Report  the  results  of  studying  the  module  characteristic  for  a  small 
SPT  and  for  a  thruster  plume.  In  particular  there  are  represented  the  description  of  test  facility  and 
system  of  thruster  performance  and  plume  parameter  measurements  as  well  as  the  results  of  theses 
measurements.  Performance  characterization  was  made  for  SPT-25  laboratory  model  having  exit 
parts  of  the  discharge  chamber  walls  made  of  the  Russian  -  borosil  type  -  ceramics  and  French  one, 
AlNBN-type,  and  of  the  USA  -  ceramics,  AX05-type. 

Plume  measurements  were  made  for  the  optimized  operation  modes  with  discharge  power 
N~100W  and  N~200W.  These  measurements  include  determination  of  the  accelerated  ion 
distribution  in  off-axis  angle  allowing  estimation  of  the  accelerated  ion  flow  divergence, 
measurements  of  the  accelerated  ions  distributions  in  energy,  plasma  potential,  electron  temperature 
and  plasma  density  number  distributions  in  off-axis  angle,  hi  the  process  of  flie  experiments,  the 
study  of  the  ion  charge  structure  in  the  SPT-25-jet  was  also  realized.  The  presence  of  one-,  two-, 
three-,  charged  Xe-  ions,  as  well  as  that  of  neutral  atoms,  was  detected. 

But  as  a  whole  the  complex  of  the  obtained  data  could  be  definitely  used  for  checking  of 
some  theoretical  interpretations  of  the  discharge  chamber  material  influence  on  thruster  operation 
and  its  characteristics. 

In  the  second  part  of  the  Report,  the  results  of  studying  some  features  of  the  physical 
processes  in  APPT,  as  well  as  an  analysis  of  the  produced  results  based  on  the  numerical  simulation 
for  increasing  their  efficiency,  are  represented.  Two-dimensional  model  of  APPT  have  been 
developed.  Modeling  shows  that  the  model  adequately  describes  formation  and  acceleration  of  a 
plasma  flow  in  coaxial  accelerating  channel.  The  model  gives  a  correct  qualitative  picture  of 
separate  stages  of  discharge  in  APPT.  The  features  of  propellant  flow  rate  determine  monotonous 
character  of  distribution  of  density,  axial  velocity  of  plasma,  and  magnetic  field  in  the  accelerating 
channel.  The  developed  numerical  model  of  APPT  is  a  tool  to  improve  thruster  parameters. 

Energy  flux  onto  propellant,  Teflon  mass  loss,  current  and  voltage,  near  propellant  bar 
electron  density  have  been  measured  in  APPT  stand.  The  energy  fraction  released  from  the 
discharge  region  upon  the  propellant  surface  in  the  APPT  is  near  10'^  of  the  energy  stored  in  the 
power  supply  source.  Main  part  of  measured  flux  is  transferred  by  particles. 

Interferometric  investigations  inside  APPT  discharge  channel,  piezoprobe  measurements  of 
dynamic  pressure  distributions  and  ultra  high-speed  photography  with  visualization  of  plasma  flows 
have  been  carried  out.  Studies  of  spatial  and  temporal  distributions  of  the  electron  density  in  the 
discharge  channel  make  possible  to  reveal  main  distinctions  of  physical  processes  in  the  high 
efficiency  APPT’s  defining  the  improved  thruster  characteristics. 

Two  APPT  models  having  the  level  of  bank  energies  20  J  and  40  J,  are  developed  and 
tested.  The  improved  APPT  characteristics  are  obtained  with  completion  of  an  optimization 
between  parameters  of  an  electric  circuit  and  the  sizes  of  an  accelerating  channel  of  the  thruster, 
when  tile  discharge  is  close  to  aperiodic. 
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1.  Introduction 

Some  results  of  studying  two  promising  diagrams  of  electrorocket  thrusters:  stationary 
plasma  -  low  power  -  thmster  (SPT)  (Task  1)  and  an  ablative  pulsed  plasma  one  (APPT)  (Task  2) 
are  given  in  this  Report.  The  studies  were  done  in  accordance  with  the  technical  calendaiy  plan  of 
the  studies  on  the  ISTC-Project  No.  2317p  for  the  period  from  April,  1, 2002  to  March,  31, 2003. 

II.  Task  1.  Study  of  stationary  plasma-low  power-thruster 


II.l.  Introduction 

As  it  is  known  /!/  the  secondary  electron  emission  plays  significant  role  in  the  physical 
processes  inside  the  SPT  accelerating  chaimel.  But  there  was  not  made  systematical  studies  of  this 
factor  impact  on  SPT  operation  and  performance  in  the  former  times.  Only  last  years  an  increase  of 
interest  to  SPT  in  the  west  countries  stimulated  some  studies  but  they  were  mainly  theoretical  ones 
/2-4/.  It  is  necessary  to  add  that  nowadays  there  are  some  plans  to  develop  the  so-called  small  SPT’s 
for  small  spacecrafts.  Because  the  less  thraster  size  the  more  difficult  to  ensure  high  its  thrust 
efficiency  it  seems  more  important  to  study  impact  of  discharge  chamber  material  for  small  SPT. 
Therefore  it  was  decided  to  study  within  the  work  program  of  the  ISTC  project  #2317p  integral  and 
plume  characteristics  of  the  SPT-25  laboratory  model  with  the  discharge  chamber  made  of  ceramics 
with  different  secondary  electron  emission. 

Particularly  it  was  necessary  to  determine: 

-  thruster  performance,  density  and  energy  distributions  of  the  accelerated  ions  at  surface  under 
control  with  the  radius  R=(0,5-0,7)m  with  center  at  the  point  of  crossing  the  thmster  axis  with 
output  plane  of  the  thmster  and  these  measurements  were  to  be  made  under  powers  ~100W  and 
~200W  for  small  SPT  with  2-3  discharge  chamber  materials; 

-  distributions  of  plasma  density,  plasma  potential,  electron  temperature  in  the  small  SPT  plume  at 
the  distances  (0,5-0,7)m  and  these  measurements  were  to  be  made  for  small  SPT  with  basic 
ceramics  under  power  ~100W; 

-  charge  stmcture  of  plasma  flux  flowing  out  from  the  thmster  in  various  directions,  neutral  fraction 
included. 

According  to  the  work  program  one  of  the  main  tasks  of  studies  within  the  ISTC  #2317 
project  was  determination  of  the  SPT-25  characteristics  for  2-3discharge  chamber  materials  having 
different  secondary  electron  emission  yield.  This  task  was  solved  during  fulfillment  of  works  and  in 
the  1**  report  161  issued  on  the  mentioned  project  there  were  represented  results  of  the  SPT-25 
characterization  in  the  cases  of  its  discharge  chamber  made  of  AIN-BN  ceramics  made  in  French 
and  BN-Si02  called  borosil  or  BGP  ceramics  produced  in  Russia. 

In  the  mean  time  at  CERN  there  were  received  data  on  these  materials  secondary  electron 
emission  yield  (SEY). 

In  the  second  half-a-year  (second  stage  of  the  studies)  the  data  on  the  SPT-25  activity 
characteristics  with  the  discharge  chamber  walls  made  of  borosil  and  of  the  ceramics,  AX05-type, 
based  on  the  BN  produced  in  the  USA  were  obtained.  Unfortunately  during  this  phase  of  tests  there 
was  broken  the  heaterless  cathode  used  in  the  U*  series  of  tests.  Therefore  there  was  realized  the 
comparative  study  of  characteristics  with  both  materials  and  new  cathode.  Nowadays  there  are  also 
available  preliminary  CERN  data  on  the  AX05  SEY.  So,  as  a  result  of  joint  efforts  coimected  with 
the  mentioned  project  fulfillment  there  was  obtained  the  U*  data  pack  for  analysis  of  three  discharge 
chamber  wall  material  and  their  SEY  influence  on  the  SPT  operation  and  characteristics.  Thus,  the 
task  mentioned  above  was  solved. 
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II.2.  Methodology  of  the  SPT-25  performance 
and  plume  characterization 

Since  early  studies  of  the  Stationary  Plasma  Thrusters  (SPT)  there  was  supposed  that  one  of 
the  key  factors  determining  the  discharge  chamber  material  (ceramics)  impact  on  thruster  operation 
is  the  secondary  electron  emission  of  ceramics  /!/.  But  as  it  was  mentioned  there  was  no  direct 
experimental  confirmation  of  this  point  as  well  as  there  was  no  systematic  study  of  the  mentioned 
impact.  Nowadays  due  to  the  wide  involvement  of  western  scientists  into  SPT  studies  there  is 
appeared  great  interest  to  the  possible  influence  of  ceramics  properties  on  plasma  dynamics  inside 
the  SPT  accelerating  chaimel  and  thruster  performance.  In  this  cormection  as  one  of  the  main  tasks 
of  study  within  the  frames  of  given  project  there  was  specified  in  a  work  program  the  study  of 
integral  characteristics  of  thrusters  having  discharge  chamber  walls  made  of  materials  with 
definitely  different  secondary  electron  emission.  These  materials  are  as  follows: 

1 .  Russian  borosil  type  ceramics  consisting  of  the  SiOi  and  BN  as  the  main  components. 

2.  AINBN  (ABN)  type  ceramics  consisting  of  BN  and  AIN  as  the  main  components. 

3.  The  ceramics,  AX05-type,  based  on  BN,  is  the  USA-produced  one. 

The  studies  were  done  at  two  stages.  At  the  first  stage,  SPT-characteristics,  using  the 
ceramics,  borosil-type  and  AIN-BN-one,  were  compared.  At  the  second  stage,  the  comparisons 
were  done  for  the  ceramics  borosil-type  and  for  tiiat  of  AX05-type.  According  to  data  obtained  at 
CERN  both  these  materials  have  significantly  higher  cross  over  energy  s'  (corresponding  to  the 
secondary  electron  emission  yield  (SEY)  1)  than  that  one  for  the  AIN-BN  ceramics  (Table  1). 

Table  1.  Cross  over  energy  s'  values  for  different  materials. 


AIN-BN 

20-40  eV 

BN-SIO2 

~80eV 

BN 

~  llOeV 

For  study  there  was  used  the  SPT-25  type  thruster  model  (Fig.  1.1.)  151  having  magnetic 
system  with  one  magnetization  coil.  Advantage  of  such  model  is  small  variation  of  the  magnetic 
field  topology  under  variation  of  magnetic  field  induction  to  optimize  operation  mode.  This  should 
simplify  analysis  of  obtained  data  presented  in  a  report.  Discharge  chamber  wall  exit  parts  (Fig.  1.2.) 
were  made  as  circular  rings  kept  to  magnetic  system  elements.  So,  it  was  possible  to  manufacture 
these  rings  of  different  materials.  As  it  was  mentioned  above  in  the  basic  option  there  was  used 
borosil  type  ceramics  delivered  by  Russian  manufacturer.  The  ABN  samples  were  delivered  by 
western  supplier  according  to  AFRL  order  and  the  second  discharge  chamber  wall  option  was 
manufactured  using  this  material.  Thruster  geometry  and  sizes  of  parts  were  identical  within  the 
manufacturing  accuracy. 

SPT-25  model  was  tested  at  test  facility  with  vacuum  chamber  of  2m  in  diameter  and  6m  in 
length  (Fig.1.3.).  This  chamber  was  equipped  by  thmstmeter  having  accuracy  of  measurements 
±3%  within  tire  range  of  the  thrust  values  5-10  mN.  Test  facility  has  also  system  supplying  gas 
(Xe)  into  the  accelerating  chaimel  through  anode  and  into  cathode  as  well  as  power  supply  sources 
for  all  thruster  circuits  and  measuring  system  (Fig.  1.4.).  The  mass  flow  rate  through  anode  was 
controlled  and  measured  with  accuracy  ±  2%  within  the  range  of  (0,5-l)mg/s  and  cathode  mass 
flow  rate  -  with  accuracy  of  ±  5%  under  its  value  0,lmg/s. 

Test  facility  was  equipped  also  by  the  system  of  the  accelerated  ion  flow  parameters 
measurement  allowing  determination  of  the  accelerated  ion  current  distribution  along  the  semicircle 
with  the  center  positioned  at  the  thruster  exit  plane  and  with  circle  plane  consisting  of  thruster  axis. 
To  realize  these  measurements  there  was  used  the  RPA  probe  mounted  on  the  boom  rotating  along 
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the  mentioned  plane  within  ±90°  relative  to  the  thruster  axis  (see  Fig.1.3.)-  This  system  allows 
estimation  of  the  flmister  plume  divergence  and  angular  distribution  of  accelerated  ions  in  energy. 

Besides  the  mentioned  main  task  there  was  specified  also  the  necessity  to  characterize 
plasma  parameter  distributions  in  a  small  SPT  plume  under  its  operation  with  discharge  power 
N(i«  lOOW  and  basic  type  of  ceramics  (borosil).  Such  characterization  was  also  made  and  obtained 
results  are  represented  below. 

To  characterize  the  tiiruster  model  performance  there  were  measured  discharge  current  and 
fimist  under  different  currents  in  magnetization  coil  for  several  fixed  mass  flow  rates  through  anode 
and  several  fixed  discharge  voltage  values.  The  range  of  mass  flow  rates  through  anode  was  varied 
within  (0,5-l,0)mg/s.  The  range  of  discharge  voltages  was  (125-250)V  and  their  values  were  varied 
with  step  AU»25V.  Thus,  it  was  possible  to  determine  possible  discharge  chamber  material  impact 
on  die  voltage-current  characteristics,  dependence  of  discharge  current  on  magnetic  field  intensity  as 
well  as  to  determine  thrust  and  thmst  efficiency  under  different  discharge  conditions. 

For  operation  modes  witii  high  enough  thrust  efficiency  and  powers  Nd=100W  and 
Nd=200W  there  was  measured  angular  distributions  of  the  accelerated  ion  current  density  and 
distributions  of  ions  in  energy  for  several  of-axis  directions. 

The  angular  accelerated  ion  current  distributions  were  determined  under  retarding  potential 
Ur=50V  relative  to  cathode  potential.  Because  plasma  potential  in  a  plume  is  (20-30)V  usage  of 
such  retarding  potential  allows  measurement  of  current  created  by  ions  with  energies  higher  than 
(20-30)eV.  It  is  necessary  to  note  that  nowadays  there  is  no  verified  methodology  for  the  accurate 
measurement  of  total  ion  current  exhausting  thruster.  Usage  of  the  mentioned  retarding  potential 
allows  characterization  of  the  angular  distribution  of  ions  in  a  plume  able  to  have  notable 
mechanical,  thermal  and  erosion  impacts  on  surfaces  crossed  by  plume. 

To  measure  the  plasma  parameter  distributions  in  a  plume  in  addition  to  RPA  there  was  used 
cylindrical  Langmuir  probe  with  axis  oriented  to  thruster  exit.  Due  to  existence  of  the  accelerated 
ions  in  a  plasma  it  is  difficult  to  have  satisfactory  interpretation  of  probe  characteristics.  Experience 
coupled  by  authors  of  this  report  shows  that  probe  characteristics  close  to  classical  one  could  be 
obtained  only  with  cylindrical  probe  oriented  as  it  was  mentioned  above.  To  verify  this  conclusion 
there  were  made  special  measurements  with  usage  of  cylindrical  probe  and  emissive  probe.  These 
measurements  gave  satisfactory  agreement  of  the  plasma  potential  values  measured  by  emissive 
probe  and  derived  from  probe  characteristics  obtained  by  cylindrical  probe. 

To  estimate  the  ion  back  flows  there  were  used  system  of  flat  probes  positioned  in  a  thruster 
exit  plane  (see  Fig.1.3.). 

Plasma  parameters  in  a  SPT-25  plume  were  determined  for  three  distances  of  probes  from 
the  thruster  exit,  namely:  along  the  semicircles  with  radiuses  R=0,3m,  R=0,4m  and  R=0,5m. 
Sensitivity  of  measuring  system  did  not  allow  obtaining  of  reliable  data  at  larger  distances. 


II.3.  Results  of  the  SPT-25  performance  characterization 

n.3.1.  SPT-25  characteristics  with  the  output  rings  made  of  the 
borosil  and  AXNBN-ceramics. 

Results  of  the  SPT-25  option  #1  (with  borosil  ceramics  exit  rings)  performance 
characterization  (Fig.  1 .5 ...  1 .28.)  show  that: 

1.  Model  has  typical  for  SPT  voltage-current  characteristics  under  great  enough  magnetic 
induction  and  there  is  saturation  of  thruster  performance  under  high  enough  magnetization 
currents.  It  is  possible  to  conclude  also  that  under  magnetization  currents  Im^(2,5-3)A  there 
is  no  further  increase  of  the  thrust  efficiency  with  increase  of  the  magnetization  current.  So, 
the  mentioned  values  could  be  considered  as  optimal  ones.  It  is  necessary  to  add  that  the 
mentioned  characteristics  were  determined  while  magnetization  current  was  reduced  from 
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high  values  till  low  ones  to  reduce  impact  of  the  thruster  model  overheating  taking  place  under 
low  magnetization  currents  while  discharge  current  is  great. 

There  is  confirmed  the  possibility  to  obtain  acceptable  total  thrust  efficiency  of  the  SPT-25 
model  with  discharge  chamber  walls  made  of  borosil  under  its  operation  with  discharge  power 
Nd*  lOOW  and  good  enough  performance  level  under  discharge  power  Nd*200W.  It  is  necessary  to 
note  that  the  thrust  efficiency  T/ta  values  represented  in  Figures  were  calculated  taking  into  account 
power  losses  for  magnetization  but  not  accounting  for  cathode  mass  flow  rate.  The  performance 
data  for  the  above  mentioned  operation  modes  parameters  are  represented  below: 


Mode  1 

Mode  2 

Mass  flow  rate  through  anode,  mg/s 

0.6 

0.7 

Mass  flow  rate  through  cathode,  mg/s 

0.1 

0.1 

Discharge  voltage,  V 

151 

Discharge  current,  A 

0.64 

0.86 

Discharge  power,  W 

96.7 

192.9 

Thrust,  mN 

6.0 

Magnetization  current,  A 

2.82 

2.93 

Voltage  at  magnetization  coil,  V 

2.49 

2.65 

Power  losses  for  magnetization,  W 

Total  power  consumption,  W 

Specific  impulse  (taking  into  account  the  cathode  mass  flow  rate) 

885 

1299 

Total  thrust  efficiency,  %  (taking  into  account  the  cathode  mass  flow 
rate  and  power  losses  for  magnetization) 

0.253 

0.324 

Plume  half  angle  for  95%  of  the  accelerated  ions  determined  by  RPA 
probe  measurements  at  R=0,5m 

78.8 

66.1 

Integral  characteristics  of  SPT-25  model  option  #3  with  discharge  chamber  walls  made  of 
AINBN  are  similar  to  that  ones  obtained  for  borosil  case  (Fig.  1.29...  1.49.)  but  under  similar 
conditions  discharge  current  is  typically  higher  (see  Fig.  1.29.,  1.36.,  1.43.)  and  thrust  is  typically 
lower.  That  is  why  thrust  efficiency  in  the  case  of  AINBN  ceramics  is  lower  (see  Fig.  1.35.,  1.42., 
1.49.).  Moreover  the  range  of  the  stable  operation  modes  is  more  narrow,  probably,  due  to  reduction 
of  the  thrust  efficiency. 

Comparing  the  voltage-current  characteristics  of  thruster  under  the  same  mass  flow  rate  one 
can  conclude  that  under  low  magnitudes  of  the  magnetization  current  the  discharge  current  in  the 
AINBN  case  is  notably  higher  and  only  under  high  magnetization  currents  the  discharge  currents  are 
comparable  for  both  cases.  So,  there  was  found  definite  difference  of  integral  characteristics  for 
different  materials.  It  is  necessary  to  add  that  under  discharge  voltages  Ud~200V  or  higher  there  is 
appeared  breakdowns  between  anode  and  internal  magnetic  pole  of  the  magnetic  system  creating 
hole  in  a  internal  discharge  chamber  exit  ring  and  with  this  hole  thmster  model  could  not  operate 
satisfactory.  After  the  1st  such  breakdown  thraster  model  was  repaired  but  after  the  2"^  breakdown 
it  was  decided  that  this  concrete  AINBN  ceramics  is  not  acceptable  for  the  SPT  discharge  chamber 
manufacturing.  Therefore  the  plume  measurements  for  this  thruster  model  option  was  not  made. 


n.3.2.  SPT-25  characteristics  with  the  output  rings  made  of  the 
borosil  and  AXOS-ceramics. 

The  main  task  of  the  given  stage  of  works  was  the  characterization  of  the  SPT-25  under  its 
operation  widi  BN  based  ceramics  of  the  AX05  type.  Unfortunately  during  new  stage  of  tests  there 
was  broken  the  heaterless  cathode  (cathode  #1)  used  in  the  1st  series  of  tests  and  there  was  no 
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possibility  to  repair  it  quickly  enough.  Therefore  it  was  decided  to  realize  new  series  of  tests  with 
the  cathode  (cathode  #2)  for  the  SPT-100  type  thruster.  Surely  this  cathode  is  oversized  for  the  SPT- 
25  and  it  is  difficult  to  operate  it  under  mass  flow  rates  through  the  cathode  c  *  0,1  mg/s. 
Therefore  there  was  used  an  additional  heating  of  this  cathode  by  power  ~25W.  Due  to  the 
mentioned  changes  the  thruster  characteristics  were  changed.  That  is  why  to  compare  these 
characteristics  for  new  pair  of  materials  there  was  made  their  determination  with  bodi  materials 
under  operation  with  new  cathode. 

As  it  was  mentioned  above  replacing  of  cathode  had  changed  the  thruster  characteristics 
significantly,  namely:  with  new  cathode  the  discharge  current  and  thrust  are  lower  than  that  ones 
with  the  heaterless  cathode  (Fig.  1.50.-1.57.)  and  finally  the  thrust  efficiency  is  also  lower 
(Fig.  1.58.-1.61.).  This  is  an  indication  that  the  SPT-25  with  new  cathode  even  with  additional 
heating  can  not  operate  effectively  enough.  Nevertheless  thruster  characteristics  are  repeatable  what 
was  checked  with  borosil  ceramics.  Therefore  it  was  decided  to  realize  new  phase  of  the 
comparative  study  of  material  impact  on  thruster  characteristics  with  new  cathode  (cathode  #2). 

The  tests  were  made  in  the  same  vacuum  chamber  of  2m  in  diameter  and  6  m  in  length,  with 
the  same  instruments  and  devices  for  the  measurement  of  electric  parameters,  thrust,  mass  flow 
rates  and  plume  parameters  as  it  was  described  in  /6/.  Additionally  there  were  measured  the 
discharge  voltage  and  current  oscillations  intensity  by  registering  of  their  (5-10)  occasional  traces  by 
digital  oscilloscope,  realizing  their  Fourier-analysis  within  the  range  of  frequencies  0-250  kHz  and 
calculating  the  mean  RMS  amplitudes  of  the  discharge  voltage  and  current  oscillations  for  50 
harmonics.  There  was  made  also  determination  of  the  dominating  oscillation  frequency.  Obtained 
oscillation  characteristics  for  the  borosil  case  confirm  that  operation  modes  with  different  cathodes 
are  different  (Fig.  1.62.-1.73.).  This  is  an  indication  of  fact  that  the  discharge  chamber  material  is 
only  one  among  many  of  factors  determining  the  SPT  operation  and  characteristics. 

There  were  made  also  the  calculation  of  the  magnetic  field  topology  (Fig.  1.74.)  using 
standard  codes  and  measurement  of  the  magnetic  field  radial  component  along  the  accelerating 
channel  mid  surface  (Fig.  1 .75.).  Obtained  results  show  that: 

-  magnetic  field  is  mainly  radial  inside  the  SPT-25  accelerating  channel  (see  Fig.  1 .74.); 

-  there  is  an  indication  of  some  magnetics  saturation  under  magnetization  currents  exceeding  ~3A 

(Fig.1.75.)  what  is  to  be  taken  into  account  under  analysis  of  the  obtained  data. 

Results  of  the  comparative  SPT-25  characteristics  determination  with  two  different  ceramics 
(Fig.  1 .76 ...  1 . 1 54.)  show  the  following: 

-  as  in  the  case  of  AIN-BN/BGP  pair  there  are  some  differences  in  voltage-current  characteristics 
under  low  magnetic  fields  in  die  case  of  borosil/AX05  pair  and  with  increase  of  magnetic  field 
these  differences  become  negligible  (Fig.  1.76. ,1.82.,  1.88.); 

-  in  the  case  of  AX05  the  discharge  current  under  low  magnetic  fields  is  typically  a  little  bit  higher 

(see  Fig.  1.76.,  1.82.,  1.88.)  and  thrust  is  a  little  bit  lower  than  in  borosil  case,  therefore  the  thrust 
efficiency  under  low  magnetic  fields  is  a  little  bit  lower  in  the  AX05  case  (see  Fig.  1.81.,  1.87., 
1.92.); 

-  dependences  of  integral  parameters  on  magnetization  current  are  qualitatively  the  same  as  for  the 
AlN-BN/borosil  pair  (see  Fig.1.93... 1.108.); 

-  the  oscillation  characteristics  (see  Fig.  1.1 09...  1.1 42.)  are  significantly  different  for  significant 
part  of  the  operation  modes  (see  Fig.  1.1 13...  1.1 15.,  1.124...  1.126.)  and  have  no  direct  link  with 
the  discharge  current  behaviour  (see  Fig.  1 .86.  and  Fig.  1.115.,  1.116); 

-  dependences  of  the  oscillation  intensity  on  he  magnetization  current  are  more  or  less  similar  for 
BGP  and  AX05  cases  (see  Fig.  1.1 27...  1.134.),  but  character  of  the  dominating  oscillation  mode 
firequency  dependences  on  this  current  are  different  at  least  under  low  voltages(see 
Fig.l. 135. ..1.142.); 
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-  the  plume  divergence  general  behaviour  is  similar  for  both  materials  (see  Fig.  1 . 1 43 ...  1 . 1 54.). 

II.3.3.  Discussion  over  the  produced  results  in  the  studies  of  an 
effect  of  various  ceramics  on  the  SPT-characteristics. 

Considering  obtained  results  it  is  possible  to  note  that  it  is  difficult  to  explain  all  of  them 
using  one  or  another  model  due  to  complexity  of  processes  in  the  SPT  discharge.  Therefore  theses 
results  could  be  considered  as  data  base  for  the  fUrther  analysis.  Taking  this  into  account  there  is 
represented  in  a  report  the  maximum  of  primary  results  because  only  complex  analysis  of  these 
results  can  give  satisfactory  their  theoretical  interpretation.  An  another  point  is  that  even  represented 
data  in  spite  of  their  great  volume  are  not  full  enough.  In  particular  there  are  no  data  on  the  local 
plasma  parameter  distributions  inside  the  accelerating  chaimel.  These  data,  especially  information 
on  the  electron  energy  distribution  along  the  accelerating  channel,  would  be  exclusively  beneficial 
for  the  analysis  of  the  discharge  chamber  material  and  its  secondary  electron  emission  influence  on 
thruster  operation  and  performance.  So,  it  seems  prospective  to  continue  study  in  order  to  get  the 
mentioned  data.  In  the  mean  time  it  seems  possible  to  make  some  conclusions,  namely: 

1 .  The  discharge  chamber  material  has  definite  impact  on  the  thruster  operation  and  performance 
(this  fact  is  well  known)  and  this  impact,  probably,  is  connected  with  the  redistribution  of  all 
local  parameters  (to  be  checked  by  the  direct  measurements)  causing  notable  change  of  thrust 
and  thrust  efficiency  even  under  great  enough  magnetic  fields  while  voltage-current 
characteristics  and  dependences  of  discharge  current  on  the  magnetization  current  are  very 
close  for  different  materials. 

The  mentioned  changes  could  be  connected  with  difference  of  the  secondary  electron 
emission  yield  but  as  it  was  mentioned  above  such  conclusion  is  to  be  checked  by  additional 
measurements.  The  main  concern  is  that  even  under  great  enough  magnetic  fields  (/„  >  2.5  A 
for  SPT-25  with  borosil  and  AIN-BN  cases)  there  is  no  notable  difference  of  the  discharge 
current  expected  for  material  with  higher  SEY  (AIN-BN)  due  to  increase  of  the  near  wall 
conductivity  but  thrust  value  in  this  case  is  evidently  less.  The  first  fact  indicate  that,  probably, 
the  ion  and  electron  currents  at  the  thruster  exit  under  comparable  operation  modes  are  close  for 
both  cases  (to  be  checked  by  measurements),  but  the  acceleration  efficiency  is  different  due  to 
difference  of  the  all  local  parameter  distributions  resulting  in  different  wall  losses,  mean  ions 
energy  etc(to  be  checked  by  measurements). 

It  is  necessary  to  add  that  for  material  with  great  secondary  emission  (cr  >  1)  under 
operational  conditions  there  is  another  possibility  to  get  higher  discharge  current  than  that  one 
coimected  with  increase  of  the  near  wall  conductivity,  namely:  there  is  the  possibility  to  obtain 
an  additional  current  through  bulk  of  insulator  as  through  conducting  material.  This  possibility 
is  confirmed  by  two  breakdowns  registered  during  the  SPT-25  operation  with  the  AIN-BN 
insulator.  These  breakdowns  appeared  during  operation  imder  high  voltages  (~200V)  and  high 
discharge  currents  (low  magnetic  fields)  and  were  providing  the  axially  directed  channel  with 
almost  circular  cross-section  and  diameter  of  ~0,5  mm,  going  in  the  internal  replacable  ring 
bulk  from  the  anode  side  to  the  cathode  side. 

Increase  of  the  discharge  current  due  to  evaporation  of  the  ceramics  material  could  be 
significant  enough  because  this  evaporation  significantly  increase  gas  medium  density  inside 
the  accelerating  channel  during  short  time.  But  in  the  prebreakdown  phase  it  is  to  be  also 
notable  to  start  significant  heating  of  the  future  breakdown  channel.  So,  probably,  an  increase 
of  the  discharge  current  under  low  magnetic  fields  for  AIN-BN  case  noted  in  /6/  was  caused  by 
this  parasite  prebreakdown  current  through  insulator  bulk.  This  current  has  occasional  character 
and  was  able  to  disturb  somewhere  the  voltage-current  characteristics  for  the  considered  case. 
This  idea  appears,  if  one  takes  into  account  that  for  many  different  comparable  operation  modes 
the  discharge  currents  are  very  close. 
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2.  If  the  discharge  current  values  for  comparable  operation  modes  and  materials  with  different 
SEY  could  be  considered  as  close  one,  this  does  not  mean  that  the  secondary  emission  and  near 
wall  conductivity  does  not  play  any  role.  Indeed,  there  were  studied  within  this  project  the 
operation  modes  with  relatively  low  discharge  voltages  and  consequently  with  not  too  high 
electron  energies.  For  these  conditions  the  mean  electron  temperatures  could  be  within  the 
range  of  10-20  eV  and  plasma  densities  near  the  walls  ~  (2  -  5)10''  1/cm^.  The  corresponding 
Debay  length  could  be  within  the  range  of  0.015-0.05  mm  that  is  comparable  with  the  scale  of 
the  surface  roughness  after  mechanical  treatment.  The  nonsmoothness  of  tiie  equipotential 
surfaces  within  the  Debay  layer  can  cause  an  additional  electron  drift  momentum  dissipation 
/?/.  This  means  that  significant  part  of  electrons  interacting  with  the  Debay  layer  can  change 
their  drift  momentum  and  give  their  input  into  the  near  wall  current  in  spite  of  the  fact  that  only 
small  part  of  electrons  was  able  to  impinge  the  wall  surface.  So,  probably,  the  part  of  electrons 
interacting  with  wall  with  appearance  of  the  secondary  emission  was  small  in  die  studied 
conditions  in  comparison  with  that  one  changing  their  drift  momentum  and  giving  impact  into 
the  near  wall  current. 

3.  Relatively  low  discharge  voltages  and  electron  temperatures  can  explain  not  so  significant 
difference  in  the  characteristics  obtained  for  the  second  studied  pair  of  materials 
(borosil/AX05).  Indeed,  due  to  low  electron  energies  there  was  not  so  much  electrons  able  to 
reach  the  level  of  energies  close  to  cross  over  energy  for  the  secondary  electron  emission. 
Probably,  due  to  this  there  was  not  appeared  significant  changes  of  thruster  characteristics  in 
spite  of  notable  difference  in  SEY.  But  considering  the  oscillation  characteristics  one  can  tell 
that  there  was  observed  definite  changes  of  thruster  operation. 

To  clarify  the  discussed  points  as  it  was  mentioned  above  it  is  necessary  to  fulfil  more 
measurements.  The  most  useful  are  to  be  the  local  plasma  parameter  measurements  inside  the 
discharge  chambers  with  walls  made  of  different  materials. 

II.4.  Results  of  the  SPT-25  plume  plasma  parameter 
measurements 

n.4.1.  Electron  temperature,  density  and  plasma  potential 
distributions  in  SPT-25  jet. 

As  it  was  mentioned  above  these  measurements  were  made  with  usage  of  RPA  and 
cylindrical  probes.  Data  obtained  by  RPA  probe  at  R=0,5m  (Fig.  1.1 55.)  gave  the  plume  half  angle 
values  represented  above.  Measurements  made  by  RPA  and  measurements  of  plasma  parameters  by 
cylindrical  probe  (Fig.l. 156.. .1.163.)  at  different  distances  gave  results  allowing  to  conclude  that: 

-  mean  energy  of  ions  is  reduced  significantly  under  off-axis  angles  higher  than  45  (see  Fig.  1.156., 
1.160.)  and  it  is  interesting  that  difference  between  discharge  voltage  and  mean  ion  energy 
maximum  is  almost  the  same  under  Ud=150V  and  Ud=225V; 

-  plasma  potential  distributions  (Fig.l. 157.,  1.161.)  has  maximum  in  vicinity  of  thruster  axis  and 
level  of  potential  (20-22V)  in  this  zone  is  close  to  that  one  obtained  for  other  thrusters; 

-  there  was  obtained  great  enough  difference  of  the  plasma  potential  values  determined  at  different 
distances  and  under  large  off-axis  angles  (seeFig.1.161.)  but  this  result  should  be  checked  in  a 
future; 

-  maximum  of  electron  temperature  is  at  level  of  ~leV  (see  Fig.l. 158.,  1.162.)  and  also  slightly 
increased  in  direction  to  the  thruster  axis; 

-  near  axis  plasma  potential  and  electron  temperature  values  are  not  changed  significantly  with 
increase  of  distance  from  the  thruster  exit; 

-  there  is  still  great  difference  (by  3-5  times)  in  plasma  density  numbers  derived  from  the  probe 
electron  current  value  at  plasma  potential  and  from  probe  ion  current  (Fig.l. 159.,  1.163.)  what 
could  be  explained  by  significant  impact  of  the  accelerated  ions  on  probe  ion  current  because  even 
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cylindrical  probe  collects  significant  number  of  these  ions  (because  thruster  is  extended  source  of 
the  accelerated  ions  moving  along  different  directions)  and  for  estimation  of  the  real  plasma 
density  number  it  seems  better  to  use  values  derived  from  the  magnitude  of  the  electron  current  at 
plasma  potential); 

-  increase  of  discharge  voltage  and  power  causes  increase  of  the  ion  current  density  in  a  near  axis 
zone  of  plume  (see  Fig.  1.155.)  and  plasma  density  number  in  this  zone  (see  Fig.1.159.  and 
Fig.  1.1 63.)  while  distributions  of  other  parameters  seem  similar  at  least  on  character  and  level  of 
magnitudes. 


n.4.2.  Plasma  charge  structure  in  the  SPT-25  jet. 

A  MX-7304  mass-spectrometer  was  used  for  measuring  the  plasma  charge  stmcture  in  the 
SPT-25  jet  in  a  given  study. 

The  equipment  disposition  in  the  measurements  is  shown  in  Fig.1.164.  The  SPT-25  was 
placed  upon  a  rotating  platform  inside  the  vacuum  chamber  which  was  evacuated  by  a  high  vacuum 
pump,  H-200,  and  the  mass-spectrometer,  MX-7304,  was  fixed  immobile  upon  a  vacuum  chamber 
flange.  The  SS-shield  protected  the  details  of  a  vacuum  lock  against  the  SPT-25  ion  beam  effect. 
The  shield  with  a  mesh  was  necessary  for  cutting  electrons  off  and  for  limiting  the  ion  beam 
(1,5mm  in  diameter).  A  potential  sufficient  for  repulsion  of  all  the  jet  electrons  was  applied  to  the 
screen  grid.  The  lens  with  a  mesh  and  the  ionizer  form  the  focusing  system  which  allows  one  to 
adjust  -  in  a  short  range  -  the  ion  beam  direction  and,  thus,  to  use  a  maximally-possible  amount  of 
the  beam  ions  for  an  analysis  with  the  mass-spectrometer.  Moreover,  a  lens  with  a  grid  could 
perform  the  role  of  an  electrostatic  ion  energy  analyzer  and,  when  necessary,  could  completely  cut 
the  ions  off,  allowing  to  pass  ttie  neutral  beam  component  only.  The  ionizer  was  used  for  an  analysis 
of  the  neutral  beam  component  and  for  the  ionization  of  a  calibrating  gas  under  the  mass- 
spectrometer  calibration  respective  to  the  mass  (mass  filter  transmission  coefficient  determination). 
MX-7304  is  a  mass-spectrometer  of  a  dynamic  type  [8],  its  principle  of  action  is  based  on  the  fact 
that  the  ions  of  a  definite  mass  (ions  of  a  substance  under  analysis),  passing  through  a  hyperbolic 
high  frequency  electric  field,  have  a  limited  amplitude  of  oscillations,  meanwhile  the  amplitudes  of 
other  ion  oscillations  rise  in  time  without  any  limits.  These  ions  enter  the  surfaces  of  electrodes, 
being  neutralized  there.  The  ions  wife  fee  limited  oscillation  amplitude  are  collected  with  fee 
collection  and  their  intensity  is  registered  -  via  an  amplifier  -  with  a  tape  recorder. 

In  the  process  of  the  experiments,  the  study  of  the  ion  structure  in  the  SPT-25  jet  was  done 
under  operation  of  fee  thruster  in  fee  power  range,  N=(l 00-200) W  at  the  Xe-consumption,  m=(0.6- 
0.7)  mg/s,  under  pressure  of  pc  =(3-4)  10‘*mm  Hg  in  the  vacuum  chamber. 

Some  fragments  of  a  typical  -  plasma-jet  -  mass-spectrogram  for  SPT-20  are  shown  in 
Fig.  1.1 65.  A  considerable  width  of  the  element  isotope  lines  is  provided  by  a  great  spread  in  fee  ion 
energies.  In  fee  treatment  of  mass  spectrograms  the  ion  current  value  with  this  or  that  charge  was 
determined  as  an  integral  over  all  the  isotopes,  i.e.  as  an  area,  S,  under  fee  intensity  distribution 
curve,y(A0>  with  fee  corresponding  coefficients,  ki  and  kj. 

■^max 

I  =  k^k2  jj(M)dM  =  k^k2S 

t^min 

Here,  the  coefficient  ki  takes  account  of  a  change  in  the  display  scale  with  respect  to  masses 
(forXe^  ki  =  1,  for  Xe^  k\  =  2);  k2  is  the  transmission  coefficient  for  fee  mass-spectrometer  probe 
with  respect  to  masses.  The  coefficient  k^  depends  on  the  mass-spectrometer  adjustment  conditions, 
and  the  precise  determination  of  this  coefficient  requires  the  mass-spectrometer  calibration  with 
respect  to  masses. 

However,  in  our  experiments,  under  mass-spectrometer  operation  conditions  with  a  low 
resolving  power,  this  coefficient  could  be  accepted  to  be  equal  k2  «1  for  Xe^  and  Xe^  /8/. 
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Here  are  some  data  on  the  plasma  charge  structure  in  the  near-axial  jet  section,  at  the  distance 
of  700mm  from  an  SPT-25  operating  under  conditions,  when  N  =  200W  (discharge  voltage,  U=225 
V),  as  an  example.  The  presence  of  one-,  two-,  three-charged  xenon  ions  was  detected,  as  well  as 
neutral  atoms,  in  the  following  proportions:  Xe^  -  63%;  Xe"^  —  13%;  Xe"^  -  3%  and  Xe°  -  21%. 

The  nature  of  some  changes  in  the  plasma  charge  structure,  dependent  on  the  SPT-25  jet 
angle,  is  shown  in  Fig.  1 . 1 66. 

The  given  dependencies  have  mainly  a  qualitative  nature,  since  the  mass-spectrometric 
useful  signal,  with  a  departure  from  the  near-axial  jet  section,  becomes  to  be  comparable  with  a 
noise  level,  accompanying  the  SPT-25  operation. 

11.5.  Conclusion 

On  base  of  the  represented  in  a  report  data  it  is  possible  to  conclude  that  the  discharge 
chamber  material  has  definite  impact  on  the  SPT  performance.  This  impact  could  be  reduced  to  the 
following: 

1.  Discharge  chamber  material  determines  the  possible  range  of  the  operation  modes.  In 
particular  it  was  obtained  that  concrete  samples  of  the  US  made  AINBN  ceramics  have 
significantly  lower  breakdown  resistively  than  borosil  type  ceramics. 

2.  Discharge  chamber  material  has  definite  impact  on  thruster  characteristics  and  performance.  In 
particular: 

-  thmster  model  with  tested  AINBN  type  ceramics  generally  ensures  larger  discharge  current  and 
lower  thrust  under  comparable  conditions  and  this  difference  is  more  significant  under  operation 
modes  with  low  magnetic  induction  inside  the  accelerating  chaimel; 

-  there  was  confirmed  the  acceptable  performance  level  of  the  SPT-25  model  with  the  borosil  type 
ceramics  (total  thrust  efficiency  ~0,25  under  specific  impulse  Isp~885s,  total  power  consumption 
N~104W  and  total  thrust  efficiency  -0,32  under  specific  impulse  Isp~1300s,  total  power 
consumption  N~201  W). 

The  results  of  measuring  the  SPT-25  plasma  jet  parameters  have  shown  that: 

1.  The  SPT-25  plume  divergence  is  higher  than  for  optimized  SPT’s  of  increased  sizes  what  is 
explained  by  usage  of  one  magnetization  coil  in  the  magnetic  system. 

2.  Plasma  parameter  distributions  in  a  SPT-25  plume  are  qualitatively  similar  to  the 
corresponding  distributions  for  larger  optimized  SPT’s  (plasma  potential  values  are  at  level  of 
-20 V,  electron  temperature  -leV).  Mean  energy  of  the  accelerated  ions  near  the  axis  are  also 
close  to  that  ones  for  the  larger  SPT’s. 

3.  In  the  SPT-25-jet,  as  well  as  in  the  SPT  of  a  greater  size,  the  presence  of  single-,  double-, 
triple-charged  xenon  ions,  as  well  as  that  of  neutral  atoms,  has  been  detected. 

As  a  result  of  works  according  to  the  work  program  there  was  obtained  the  data  pack 
allowing  together  with  results  of  the  secondary  electron  emission  determination  to  realize  some 
analysis  of  this  emission  impact  on  the  SPT  operation  and  characteristics.  It  is  shown  also  the 
necessity  of  the  further  studies  and  one  of  the  prospective  ways  to  clarify  the  appeared  questions  is 
to  realize  the  local  plasma  parameter  measurements  inside  discharge  chambers  with  wall  made  of 
materials  with  different  SEY.  Thus,  results  of  the  fulfilled  works  give  new  data  on  the  discharge 
chamber  material  impact  on  the  SPT  operation  and  performance. 
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Fig.1.1.  SPT-25  experimental  unit 
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Fig.1.3.  Schematic  of  the  vacuum  chamber. 
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Fig.1.5.  Discharge  current  versus  the  magnetization  current 
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SPT-25  #1,  ma=0.5mg/s 

N,W 


◄ - 

— ^Ud  =  150V  -*-Lld  =  175V  — *-Ud  =  200V  — ^Ud=225V  -*-Ud=250V 

Fig.1.6.  Discharge  power  versus  the  magnetization  current 


25 


ISTC 


FINAL  TECHNICAL  REPORT 


Partner  Pro  ject  i^231 7v 
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Fig.1.7.  Thrust  versus  the  magnetization  current 
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SPT-25  #1,ma=0.5  mg/s 

11  ta 


◄ - 

— ^Ud  =  150V  Ud  =  175V  ^►-Ud  =  200V  — ^Ud  =  225V  -H»-Ud  =  250V 

Fig.1.8.  Thrust  efficiency  versus  the  magnetization  current 
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SPT-25  #1,  ma=0.6mg/s 


-^Ud  =  125V  -*-Ud  =  150V  -<^Ud  =  175V 

-^Ud  =  200V  -^Ud  =  225V  -*-Ud  =  250V 


Fig.1.9.  Discharge  current  versus  the  magnetization  current 
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SPT-25  #1,ma=0.6mg/s 
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-^Ud  =  200V  -«-Ud  =  225V  -Hii-Ud  =  250V 

Fig.1.10.  Discharge  power  versus  the  magnetization  current 
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SPT-25  #1,  ma=0.6  mg/s 

T,  mN 


◄ - 

— ^Ud  =  125V  -*-Ud  =  150V  -*-Ud  =  175V 

-^Ud  =  200V  -^Ud  =  225V  -■-Ud  =  250V 


Fig.1.11.  Thrust  versus  the  magnetization  current 
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Fig.1.12.  Thrust  efficiency  versus  the  magnetization  current 
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SPT-25  #1,  ma=0.7mg/s 

N,  W 


◄ - 

— ^Ud=125V  -^Ucl  =  150V  -*-Ud  =  175V 

— ^Ud  =  200V  -«-Ud  =  225V  -*^-Ud  =  250V 

Fig.1.14.  Discharge  power  versus  the  magnetization  current 
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-^Ud  =  200V  -«-Ud  =  225V  -«-Ud  =  250V 

Fig.1.16.  Thrust  efficiency  versus  the  magnetization  current 
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SPT-25  #1,  ma=0.8mg/s 

N,  W 


◄ - 

— ^Ud=125V  -*-Ud  =  150V  -<^Ud  =  175V 

-*-Ud  =  200V  -H-Ud  =  225V  -«i-Ud  =  250V 


Fig.1.18.  Discharge  power  versus  the  magnetization  current 
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Fig.1.20. 
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Discharge  current  versus  the  magnetization  current 
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SPT-25  #1,  ma=0,9mg/s 
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— ^Ud  =  125V  -^Ud  =  150V  -*-Ud  =  175V 

-^Ud  =  200V  -*«-Ud  =  225V 


Fig.1.22.  Discharge  power  versus  the  magnetization  current 
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Fig.1.23.  Thrust  versus  the  magnetization  current 
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SPT-25  #1,  ma=1.0mg/s 
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Fig.1.26.  Discharge  power  versus  the  magnetization  current 
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Fig.1.28.  Thrust  efficiency  versus  the  magnetization  current 
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SPT-25,  nna=0.6  mg/s 

Id,  A 


-♦-#1,  Im  =  1.0  A 
Im  =  1.5  A 
--^#1,  Im  =  2.5  A 
—•-#1,  Im  =  4.0  A 

Fig.1.29.  Voltage-current  characteristics  of  the  SPT-25  with  discharge 
chamber  exit  rings  made  of  borosil  (#1)  and  ABN  (#3) 


-^#3,  Im  =  1.0  A 
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-o-#3,  Im  =  4.0A 
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Fig.1.30.  Thrust  versus  the  discharge  voltage 
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SPT-25,  ma=0.6  mg/s 
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Thrust  versus  the  discharge  voltage 
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Thrust  efficiency  versus  the  discharge  voltage 
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Fig.1.36.  Discharge  current  versus  the  discharge  voltage 
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SPT-25,  nna=0.7  mg/s 
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-•-#1,  Im  =  1.5A  -B-#3,  lm  =  1.5A 

Fig.1.38.  Thrust  versus  the  discharge  voltage 
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Fig.1.39. 
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Fig.1.40.  Thrust  versus  the  discharge  voltage 
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Fig.1.41.  Discharge  power  versus  the  discharge  voltage 
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Fig.1.42.  Thrust  efficiency  versus  the  discharge  voltage 
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Fig.  1.43.  Discharge  current  versus  the  discharge  voltage 
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Fig,l,46.  Thrust  versus  the  discharge  voltage 
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Fig.1.48.  Discharge  power  versus  the  discharge  voltage 
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Fig.1.50.  Voltage-current  characteristics  with  cathode  #1  and  cathode  #2, 
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Fig.1.51.  Voltage-current  characteristics  with  cathode  #1  and  cathode  #2. 
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Fig.1.53.  Voltage-current  characteristics  with  cathode  #1  and  cathode  #2. 
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Fig.1.54,  Thrust  versus  the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.55.  Thrust  versus  the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.56,  Thrust  versus  the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.57.  Thrust  versus  the  discharge  voltage  with  cathode  #1  and  cathode  #2. 


76 


ISTC 


FINAL  TECHNICAL  REPORT 


Partner  Pro  ject  U23 1 7p 


Tit 


SPT-25,  lm=1.0A 


280  Ud,V 


-BGP#1,  m=0.5  mg/s 
“BGP#1,  m=0.6  mg/s 
■BGP#1,m=0.7mg/s 
-BGP#1,  m=0.8  mg/s 


-0-BGP#2,  m=0.5  mg/s 
HBh~  BGP#2,  m=0.6  mg/s 
-A-BGP#2,  m=0.7  mg/s 
--e—  BGP#2,  m=0.8  mg/s 


Fig.  1.58.  Thrust  efficiency  versus  the  discharge  voltage 

with  cathode  #1  and  cathode  #2. 
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Fig.  1.59.  Thrust  efficiency  versus  the  discharge  voltage 

with  cathode  #1  and  cathode  #2. 
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Fig.  1.60.  Thrust  efficiency  versus  the  discharge  voltage 

with  cathode  #1  and  cathode  #2. 
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Fig.1.61.  Thrust  efficiency  versus  the  discharge  voltage 
with  cathode  #1  and  cathode  #2. 
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Fig.1.64.  The  discharge  voltage  oscillations  RMS  amplitude  versus 
the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.65.  The  discharge  current  oscillations  RMS  amplitude  versus 
the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.66.  The  discharge  voltage  oscillations  RMS  amplitude  versus 
the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.67.  The  discharge  current  oscillations  RMS  amplitude  versus 
the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.69. 


The  discharge  current  oscillations  RMS  amplitude  versus 
the  discharge  voltage  with  cathode  #1  and  cathode  #2. 
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Fig.1.70.  The  frequency  of  the  dominating  discharge  voltage  (FO-U)  and 

discharge  current  (FO-I)  oscillation  mode  versus  the  discharge  voltage 
with  cathode  #1  and  cathode  #2. 


89 


ISTC 


FINAL  TECHNICAL  REPORT 


Partner  Pro  ject  ^231 7p 


SPT-25,  m=0.6  mg/s,  Im=2.7  A . 


FO,  Hz 


-■-FO-U,  BGP#1 
-♦-FO-I,  BGP#1 


FO  -  U,  BGP#2 
FO  - 1.  BGP#2 


Fig.1.71.  The  frequency  of  the  dominating  discharge  voltage  (FO-U)  and 

discharge  current  (FO-I)  oscillation  mode  versus  the  discharge  voltage 
with  cathode  #1  aud  cathode  #2. 
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Fig.1.72.  The  frequency  of  the  dominating  discharge  voltage  (FO-U)  and 

discharge  current  (FO-I)  oscillation  mode  versus  the  discharge  voltage 
with  cathode  #1  and  cathode  #2. 


91 


_ Partner  Project  #2317i 

FINAL  TECHNICAL  REPORT 


SPT-25,  in=0.7  mg/s,  Im=2.7  A . 

FO,  Hz 


100  120  140  160  180  200  220  240  260  Ud,  V 


FO  -  U,  BGP#1  -S-  FO  -  U,  BGP#2 

-♦-FO-I,  BGP#1  -0-FO-l.  BGP#2 


Fig.1.73.  The  frequency  of  the  dominating  discharge  voltage  (FO-U)  and 

discharge  current  (FO-I)  oscillation  mode  versus  the  discharge  voltage 
with  cathode  #1  and  cathode  #2. 
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Fig,1.76.  Comparative  voltage-current  characteristics  with  cathode  #2. 
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Fig.1.77.  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.78.  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.81.  Thrust  efficiency  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.82.  Comparative  voltage-current  characteristics  with  cathode  #2, 
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Fig.1.83.  Thrust  versus  the  discharge  voltage  vdth  cathode  #2. 
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Fig.1.84,  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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ig.1.85.  Thrust  versus  the  discharge  voltage  with  cathode  #2, 
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ig.1.86.  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.87.  Thrust  efficiency  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.88.  Comparative  voltage-current  characteristics  with  cathode  #2 
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Fig.1.89.  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.91,  Thrust  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.1.92.  Thrust  efficiency  versus  the  discharge  voltage  with  cathode  #2. 
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Fig.  1.95.  Discharge  current  versus  the  magnetization  current. 
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Fig.  1.97.  Discharge  current  versus  the  magnetization  current. 
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1.100.  Thrust  efficiency  versus  the  magnetization  current. 
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Fig.1.101.  Discharge  current  versus  the  magnetization  current. 
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Fig.1.102.  Discharge  current  versus  the  magnetization  current. 
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Fig.1.104.  Discharge  current  versus  the  magnetization  current. 
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Fig.1.105.  Discharge  current  versus  the  magnetization  current. 
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Fig.U06.  Discharge  current  versus  the  magnetization  current. 
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Fig.1.107.  Thrust  versus  the  magnetization  current. 
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Fig.1.108.  Thrust  efficiency  versus  the  magnetization  current. 
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Fig.1.109.  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.  1.1 10.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.1 1 1 .  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.112.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.113.  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.114.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.115.  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.116.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.117.  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.118.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.U19.  The  discharge  voltage  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.120.  The  discharge  current  oscillations  RMS  amplitude 

versus  the  discharge  voltage. 
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Fig.1.121.  Dominating  discharge  voltage  (FO-U)  and  discharge  current  (FO-I) 
oscillations  mode  frequency  versus  the  discharge  voltage. 
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Fig.1.122.  Dominating  discharge  voltage  (FO-U)  and  discharge  current  (FO-I) 
oscillations  mode  frequency  versus  the  discharge  voltage. 
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Fig.1.123.  Dominating  discharge  voltage  (FO-U)  and  discharge  current  (FO-I) 
oscillations  mode  frequency  versus  the  discharge  voltage. 
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Fig.1.124.  Dominating  discharge  voltage  (FO-U)  and  discharge  current  (FO-I) 
oscillations  mode  frequency  versus  the  discharge  voltage. 
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Fig.L127.  The  discharge  voltage  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.128.  The  discharge  voltage  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.129.  The  discharge  current  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.130.  The  discharge  current  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.131.  The  discharge  voltage  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.132.  The  discharge  voltage  oscillations  RMS  amplitude 
versus  the  magnetization  current. 


ISTC 


FINAL  TECHNICAL  REPORT 


Partner  Project  #2317d 


SPT-25  B6P#2,  m=0.7  mg/s. 


-Id,  A 


-♦-Ud=150V  -»-Ud=175V  -A-Ud=200V  -•-Ud=225V  -*-Ud=250V 


Fig.1.133.  The  discharge  current  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.134.  The  discharge  current  oscillations  RMS  amplitude 
versus  the  magnetization  current. 
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Fig.1.135.  The  dominating  discharge  voltage  oscillation  mode  frequency 

versus  the  magnetization  current. 
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Fig.1.136.  The  dominating  discharge  voltage  oscillation  mode  frequency 

versus  the  magnetization  current. 


Im,  A 


154 


JSTC 


FINAL  TECHNICAL  REPORT 


Partner  Protect  #231  7d 


SPT-25  BGP#2,  m=0.6  mg/s. 

FO,  Hz 


-•-Ud=150V  -•-Ud=175V  -A-Ud=200V  -♦-Ud=225V  ^l{-Ud=250V 


Fig.1.137.  The  dominating  discharge  current  oscillation  mode  frequency 

versus  the  magnetization  current. 
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Fig.1.138.  The  dominating  discharge  current  oscillation  mode  frequency 

versus  the  magnetization  current. 
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Fig.1.139.  The  dominating  discharge  voltage  oscillation  mode  frequency 

versus  the  magnetization  current. 
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Fig.1.140.  The  dominating  discharge  voltage  oscillation  mode  frequency 

versus  the  magnetization  current. 
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SPT-25  BGP#2,  m=0.7  mg/s. 

FO,  Hz 


-♦-Ud=150V  -■-Ud=175V  -A-Ud=200V  -•-Ud=225V  -*-Ud=250V 


Fig.1.141.  The  dominating  discharge  current  oscillation  mode  frequency 

versus  the  magnetization  current. 
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SPT-25  BN#2,  m=0.7  mg/s. 

FO,  Hz 


-♦-Ud=150V  -•-Ud=175V  -Tlr-Ud=200V  -•-Ud=225V  -*-Ud=250V 


Fig.1.142.  The  dominating  discharge  current  oscillation  mode  frequency 

versus  the  magnetization  current. 
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Po.95  degrees 


SPT-25,  m=0.6  mg/s,  lm=1.5  A 


100  120  140  160  180  200  220  240  260  Ud,  V 


-•-BGP#2  -■-BN#2 


Fig.1.144.  Plume  half  angle  versus  the  discharge  voltage. 
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degrees 


SPT-25,  m=0.7  mg/s,  lm=1.5  A 


7 


100  120  140  160  180  200  220  240  260  Ud,  V 

-♦-BGP#2  -■-BN#2 


Fig.1.145.  Plume  half  angle  versus  the  discharge  voltage. 
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Po.95  degrees 


SPT-25,  m=0.8  mg/s,  lm=1.5  A 


100  120  140  160  180  200  220  240  260  Ud,  V 


-♦-BGP#2 


Hil-BN#2 


Fig.1.146.  Plume  half  angle  versus  the  discharge  voltage. 
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Po.95  degrees 


SPT-25,  m=0.5  mg/s,  lm=2.7  A 


-♦-BGP#2  -■-BN#2 


Fig.1.147.  Plume  half  angle  versus  the  discharge  voltage. 
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Po.95  degrees 


SPT-25,  m=0.6  mg/s,  lm=2.7  A 


100  120  140  160  180  200  220  240  260  Ud,  V 


-•-BGP#2  -«-BN#2 


Fig.1.148.  Plume  half  angle  versus  the  discharge  voltage. 
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Po.95  degrees 


SPT-25,  m=0.7  mg/s,  lm=2.7  A 


100  120  140  160  180  200  220  240  260  Ud,V 


-♦-BGP#2 


-•-BN#2 


Fig.1.149.  Plume  half  angle  versus  the  discharge  voltage. 
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Fig.1.150.  Plume  half  angle  versus  the  discharge  voltage. 
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Po.95  degrees 


SPT-25,  m=0.8  mg/s,  lm=4.0  A 


100  120  140  160  180  200  220  240  260  Ud,  V 


■BGP#2 


-»-BN#2 


Fig.1.154.  Plume  half  angle  versus  the  discharge  voltage. 
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SPT-25  #1, 1=0.5  m 

Ji,  mA/cm* 


-♦-N=100W  -  N=100W  -^N=200W  - N=200  W 

Fig.1.155.  The  accelerated  ion  current  distributions  in  off-axis  angle 
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SPT-25  #1.  N=100  W,  Ud=150  V,  ma=0.6  mg/s. 

Ei,  eV 


— Lpr=0.3  m  —a—  Lpr=0.4  m  — Lpr=0.5  m 
Fig.1.156.  Mean  energy  of  ions  versns  the  off-axis  angle 
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SPT-25  #1.  N=100  W,  Ud=150  V.  ma=0.6  mg/s. 

UO,  Upl,  V 


— UO,  Lpr=0.3  m  — ^  Upl,  Lpr=0.3  m 

— ^  UO,  Lpr=0.4  m  —A—  Upl,  Lpr=0.4  m 

— UO,  Lpr=0.5  m  — e—  Upl,  Lpr=0.5  m 


Fig.1.157.  Probe  floating  potential  Uq  and  plasma  potential  Upi  values 

versus  the  off-axis  angle 
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SPT-25  #1,  N=100  W,  Ucl=150  V,  ma=0.6  mg/s. 


Te.eV 


— Lpr=0.3  m  —*r-  Lpr=0.4  m  — Lpr=0.5  m 


Fig.1.158.  Electron  temperature  versus  the  off-axis  angle 
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SPT-25  #1,  N=100  W,  Ud=150  V,  ma=0.6  mg/s. 


Fig.1.159.  Plasma  density  numbers  derived  from  ion  current  value  (Ni) 
and  electron  current  value  (Ne)  versus  the  off-axis  angle 
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SPT-25  #1,  N=200  W,  Ud=225  V,  ma=0.7  mg/s. 

Ei,  eV 


-90  -60  -30  0  30  60  90 

— Lpr=0.3  m  — Lpr=0.4  m  — Lpr=0.5  m 
Fig.1.160.  Mean  energy  of  ions  versns  the  off-axis  angle 


178 


SPT-25  #1,  N=200  W.  Ud=225  V,  ma=0.7  mg/s. 

UO,  Upl,  V 


— UO,  Lpr=0.3  m  — o—  Upl,  Lpr=0.3  m 

—A—  UO,  Lpr=0.4  m  -t-  Upl,  Lpr=0.4  m 

— UO,  Lpr=0.5  m  — Upl,  Lpr=0.5  m 


Fig.1.161.  Probe  floating  potential  Uo  and  plasma  potential  Upi 

values  versus  the  off-axis  angle 
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SPT-25  #1,  N=200  W.  Ud=225  V,  ma=0.7  mg/s. 


Te.eV 


— Lpr=0.3m  — *^Lpr=0.4m  — Lpr=0.5m 


Fig.1.162.  Electron  temperature  versus  the  off-axis  angle 
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SPT-25  #1,  N=200  W,  Ud=225  V,  ma=0.7  mg/s. 


Ni,  Ne,  1/cm^ 
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Fig.1.164.  Circuit-diagram  for  measuring  the  mass  (charge) 
composition  of  the  SPT-25-plasma  jet  with 
a  mass-spectrometer,  MX  7304. 
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Xe° 

63% 

BUM 

3% 

21% 

Fig.1.165.  Mass-spectrogram  of  SPT-25-plasma  jet 
(N  =  200W,  U  =  225V,  Pch=  4  10  ®  mm  Hg). 
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Fig.1.166.  Distribution  of  ion  components,  Xe^  and  Xe^, 

dependent  on  the  angle  in  the  SPT-25  jet. 
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III  Task  2.  Ablative  Pulsed  Plasma  Thruster  Study 

III.l  Introduction 


Ablative  Pulsed  Plasma  Thrusters  (APPT)  were  the  first  application  of  electric  propulsion  in 
space  almost  40  years  ago  The  basic  operation  of  the  APPT  consists  of  repeated  discharge 
pulses  across  a  solid  propellant  surface.  Useful  thrust  is  mainly  produced  by  the 
electromagnetic  acceleration  of  the  ablated  mass  that  has  been  ionized.  The  similar  technique 
automatically  provided  the  matching  of  a  propellant  feed  with  the  APPT  parameters  and 
allows  one  to  produce  relatively  effective  plasma  acceleration.  The  applications  ranged  from 
control  propulsion  for  larger  satellites  to  primary  propulsion  for  small  satellites.  Currently, 
APPTs  are  considered  as  an  attractive  propulsion  option  for  stationkeeping  and  drag  makeup 
purposes  of  mass  and  power  limited  satellites. 

However,  the  potential  use  of  APPTs  for  formation-flying  applications  require  longer-life 
components  and  lower  system  mass  to  meet  mission  requirements.  Moreover,  application  of 
APPTs  to  orbit  raising  maneuvers  of  power-limited  spacecraft  require  higher  performance 
than  has  been  previously  demonstrated.  External  circuit  resistance,  mass  loss  at  low  speed, 
non-matched  propellant  flow  rate  multiply  to  provide  a  total  thruster  efficiency  of  less  than 
10%.  By  applying  each  of  the  inefficiencies  in  turn,  it  should  be  possible  to  improve  the 
performance  of  the  PPT  substantially 

In  the  90’s,  the  low  thrust  capabilities  of  APPT  have  been  reopened,  thanks  to  the  tendency  to 
decrease  satellite  mass.  As  the  same  time  the  experiments  on  improving  of  APPT 
performance  have  given  noticeable  results.  Russian  APPTs  demonstrated  efficiency  on  the 
level  20  -30%  for  bank  energy  range  60  -  150  J  for  the  last  few  years.  APPT  with  such 
high  thrust  efficiency  make  possible  solving  of  many  control  problems  for  small  spacecrafts 
(SSC).  This  opens  up  new  fields  of  their  use  as  propulsion  systems  for  SSC  orbit  control 
Modeling  shows  that  the  APPTs  can  enhance  their  efficiency,  when  one  manages  to 
control  the  propellant  flow  rate,  because  the  discharge  parameters  are  strongly  dependent  on 
boundary  conditions  in  the  inlet  of  accelerating  channel.  Therefore  the  understanding  of 
propellant  ablation  and  ionization  is  very  important  for  adequate  numerical  simulation  of  a 
discharge  in  APPT  acceleration  channel  and  thrusters  development  Modeling  of  APPT 
using  experimentally  obtained  dependencies  of  propellant  flow  rate  presents  useful  approach 
for  numerical  modeling  of  APPT.  However,  it  is  necessary  to  have  adequate  understanding  of 
propellant  ablation  and  ionization  at  the  inlet  of  APPT  to  control  this  for  APPT  discharge 
optimization^^l 

According  Work  Plan  (Task  2),  this  section  presents: 

Development  of  numerical  APPT  model. 

Study  of  heat  flux  transfer  to  propellant  and  ablation  losses. 

Study  of  plasma  parameters  and  discharge  optimization. 

Development  and  tests  of  the  APPT  model. 

We  developed  the  2D  MHD  model  of  plasma  acceleration  in  coaxial  APPT  using 
experimental  data.  Main  attention  is  devoted  to  conditions  at  the  inlet  of  accelerating  channel, 
which  dramatically  determine  plasma  dynamics  in  accelerating  channel.  Boundary  conditions 
for  propellant  flow  rate  are  based  on  the  of  energy  flux  propagated  from  a  discharge  to 
propellant.  Propellant  flow  rate  will  be  accounted  from  kinetics  of  thermal  degradation  of 
Teflon  with  use  of  measured  energy  flux  onto  propellant.  Data  obtained  will  be  analyzed, 
approximated  to  be  serving  as  boundary  condition  for  APPT. 
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Experimental  efforts  were  intended  to  the  study  of  energy  fluxes  transferred  from  plasma  onto 
propellant  surface.  Experiments  were  produced  on  APPT  stand  at  Kurchatov  Institute.  The 
method  of  energy  flux  measurement  and  obtained  results  are  presented.  Energy  flux 
measurements  together  with  Teflon  mass  loss  and  current  and  voltage  measurements  have 
been  carried  out.  Experiments  to  split  particle  flux  and  radiation  flux  are  described  and 
discussed. 

Improvements  in  APPT  performances  have  been  done  due  to  finding  the  optimum  relations 
between  fire  electric  circuit  parameters  and  the  discharge  channel  dimensions,  realizing  the 
quasi-aperiodic  discharge  in  the  thruster  and  optimizing  the  current  distribution  in  the 
discharge  chaimel 

For  further  APPT  efficiency  increasing  as  well  as  for  retaining  the  achieved  efficiency  under 
APPT  flight  prototypes  manufacturing  it  is  necessary  to  understand  more  clearly  the  physical 
processes  in  thruster  discharge,  chaimel.  For  this  purpose  the  interferometric  investigations 
inside  APPT  discharge  chaimel,  piezoprobe  measurements  of  dynamic  pressure  distributions 
and  ultra  high-speed  photography  with  visualization  of  plasma  flow  were  undertaken. 
Besides,  two  APPT  mock-ups  operating  at  energy  levels  of  20  J  and  40  J  were  developed  and 
tested. 

The  objective  of  this  work  is  to  analyze  the  distinctive  features  of  physical  processes  in  APPT 
to  study  these  through  the  use  of  different  diagnostic  methods,  tests  and  numerical  modeling. 
This  makes  possible  to  develop  next  generation  APPT,  with  increased  thrust  characteristics. 


III.2  Semi-Phenomenological  Numerical  Model  of  Plasma  Acceleration  in 
the  APPT  Thruster. 

III.2.1  Problem  statement 

A  rather  complete  model  of  a  plasma  flow  in  APPT  should  take  account  of  the  vapour 
molecule  dissociation  and  neutral  particle  ionization;  the  model  should  represent  the  dynamics 
of  molecules,  atoms  ions  and  electrons.  However,  even  some  attempts  to  represent  the  gas 
ionization  only  in  a  classical  approach  of  the  stationary  flow  encounter  essential  difficulties, 
since  the  gas  ionization  process  in  a  thin  layer  is  non  -equilibrium  one  and  the  temperatures  of 
ions  and  electrons  can  essentially  differ  from  each  other.  Moreover,  one  should  take  account 
of  the  pulsed  insulator  evaporating  kinetics.  In  this  situation,  it  often  turns  out  to  be  expedient 
to  use  a  power  approach  to  the  description  of  an  ionization  process  in  plasma,  using  the  idea 
of  an  electron  cost,  e 


s  =  I+3/2Te+R, 

where  I  is  the  ionization  potential,  Te  is  the  electron  temperature,  R  is  radiation  energy 
fraction  per  one  act  of  ionization.  The  introduction  of  8  into  consideration  essentially 
simplifies  the  gas  ionization  process  consideration  in  APPT  On  the  other  hand,  the 
attempts  of  strict  and  sequential  ionization  process  consideration  are  often  not  confirmed  by 
the  experiments,  since  the  development  of  ionization  and  radiation  are  mainly  determined  by 
the  presence  of  impurities  in  the  plasma.  Rather  small  plasma  temperature  changes  in 
acceleration  channel  are  a  specific  feature  of  the  electric  discharge  in  APPT.  This  fact  is 
related  with  a  strong  plasma  radiation  intensity  growth  with  an  increase  in  temperature.  Weak 
temperature  changes  in  the  acceleration  channel  allow  one  to  use  an  isothermal  model  of  the 
acceleration. 

Then,  the  plasma  dynamics  in  APPT  can  be  represented  by  the  set  of  magnetic 
hydrodynamics  equation,  which  has  the  following  form  in  CGS  -  units: 
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+  divpv  =0 


av  VP  1  r  n 

— +  (vV)v  =  - - +  —  rotH.H 

at  ^  p  47tp^  J 

— = - AH+  ro^V,  Hjf  — —rot^ - — 

dt  Ana  Ane  p 


P  = 


2pkT 

M, 


T = Tj  =  Tj  =  const 


(1) 

(2) 

(3) 

(4) 

(5) 


The  generally-  adopted  designations  of  the  values  are  used  in  Equations.  (1-5):  p  is  the 
plasma  density,  v  is  the  velocity,  P  is  the  pressure,  H  is  the  magnetic  field,  T  is  the 
temj)erature,  ct  is  the  conductivity,  Mi  is  the  ion  mass,  c  is  the  velocity  of  light,  e  is  the 
electron  charge.  In  the  equation  of  motion  we  neglect  the  viscosity  forces;  in  the  equation  of 
magnetic  field  induction,  the  electron  pressure  gradient.  The  last  term  in  the  equation  of 
induction  (3)  takes  account  of  the  Hall  effect  in  plasma. 


III.2.2  Two-dimensional  non-steady  model 


Experiments  have  shown  that  the  plasma  flow  in  coaxial  APPT  are  essentially  two- 
dimensional  ones;  there  are  current  loops  in  the  accelerating  channel,  flow  focusing  towards 
the  APPT-axis  is  observed  etc.  Therefore  a  two  dimensional  model  is  expedient  to  be  used  for 
an  adequate  representation  of  a  plasma  flow  geometry  and  focusing  in  the  instrument  under 
consideration.  APPT  is  schematically  represented  in  Figure  III.l. 


Figure  III.l.  APPT  and  the  electric  circuit  loop;  1-  internal  electrode,  2  -  propellant  bar,  3  - 
external  electrode,  4  -plasma. 


I1I.2.3  Equations 

Equations,  representing  the  plasma  acceleration  within  the  framework  of  a  two- 
dimensional  model,  can  be  obtained  from  Eqs.  (1-5)  and  from  the  relationships  for  the  initial 
and  boundary  conditions  considered  above. 

In  this  case,  as  in  one-dimensional  model,  in  Ohm's  law  for  the  plasma,  the  electron 
pressure  gradient,  the  thermal  force,  and  the  Hall  term  are  dropped.  Moreover,  it  assumed  that 
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the  degree  of  plasma  ionization  is  c  =  1,  and  T=Te=Ti.  The  plasma  conductivity  is  calculated 
from  the  relation  ct  =  1.5  lO'^  (the  Coulomb  logarithm  =6).  It  is  also  assumed  that  the 
plasma  temperature  is  constant  at  3  eV.  This  last  assumption  distorts  the  real  distribution  of 
values  in  the  acceleration  channel,  especially  at  the  leading  edge  of  the  stream  and  at  the 
dielectric  surface. 


II1.2.4  Dimensionless  equations 

The  choice  of  main  units  is: 
tQ=l/c-jL^ ,  -  characteristic  time  of  electric  circuit; 

Ro=Ri-  radius  of  central  electrode  of  accelerating  channel; 

Vo  is  tile  initial  voltage  on  the  capacitor  bank. 

Other  units  are: 

Velocity  -  vo=Ro/to,  magnetic  field  is  Ho=cVoto/Ro^,  plasma  density  is  po=Ho^/47tvo^,  energy 
unit  is  Eo=CoVoV2,  current  is  cRoHo/2,  resistance  unit  -  rgo=2vo/c^. 

Taking  the  assumption  made  above  into  account,  we  write  a  system  of  dimensionless 
equation  with  the  same  notation  for  the  quantities  used  in  the  problem: 

Then  dimensionless  equations  can  be  written  as: 
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at  RaR  dz 
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,  Morozov  parameter 


I1I.2.5  Boundary  and  initial  conditions 

in.2.5.1  Inlet  plasma  density 

Propellant  flow  rate  is  determined  by  dissipation  of  energy  transferred  from  electrical 
discharge  and  can  be  obtained  from  approach  of  Teflon  ablation  with  accounting  of  kinetics 
(see  Appendix  1). 
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In  considered  problem  statement,  boundary  condition  at  the  inlet 
is  referred  to  the  outer  part  of  ionization  layer  where  ionization 
degree  is  near  1 .  We  neglect  the  processes  in  thin  layer  where 
plasma  is  ionized.  Only  accoimt  of  energetic  losses  for 
ionization  will  take  place.  Propellant  flow  rate  will  be  accounted 
according  the  procedure  shown  in  Figure  III.2. 


Figure  III.2.  Inlet  propellant  flow  rate  obtaining  procedure 


Boundary  condition  for  propellant  (Teflon)  flow  rate  is  based  on 
energy  flux  measurements  and  approximation.  Data  obtained 
will  be  used  for  solving  a  problem  of  Teflon  destruction  and 
ablation  due  to  energy  fltix  propagated  from  electrical  discharge. 
The  version  of  the  device  for  registration  of  energy  flow  falling  on  APPT  propellant  has  been 
developed.  The  basis  of  the  registration  is  the  thermistor  (1.6 *3. 2  mm^,  weight  9  mg), 
assembled  in  a  bridge  scheme.  The  estimations  show,  that  the  sensitivity  of  the  sensor  is 
sufficient  for  measurement  of  energy  flows  on  propellant  surface. 


The  boundary  condition  for  magnetic  field  are  determined  by  the  equations  of  electric  circuit 
of  the  APPT: 

Electric  circuit  equation  can  be  written  as: 


V  -ln(A)[v^/f(l,0,0  +  ^2^i/(l,z,0|,=o]-^l^^(lA0-^(lA0'',  =0 

^  =  (7-8) 

at 


Al=  L/2Ro-  ratio  of  initial  inductance  to  diameter  of  central  electrode,  A,=R2/Ri. 

At  the  initial  stage  of  the  discharge,  the  gap  between  the  APPT-electrodes  is  filled  with  some 
amount  of  plasma  due  to  an  initiating  device.  As  a  rule,  the  relative  amount  of  igniting 
plasma,  at  various  discharge  ignition  techniques,  is  no  more,  the  a  few  percent  of  a  total 
plasma  mass,  taking  part  in  the  discharge.  Therefore  one  can  assume  that  the  plasma  occupies 
a  cylindrical  layer,  1  -thick,  at  t=0: 

0<z<z*,  Ri<R<R2, 

and  it  is  characterized  by  the  following  parameters 

Vz=Vo,  V,  =0,  H  =  0,  T=To,  p  =  p* 


p*  -  is  initial  plasma  density  in  thin  cylindrical  layer  at  t=0.  Value  of  p*  provide  initial 
plasma  mass  at  the  level  of  (1-10)%  of  total  evaporated  mass  per  discharge. 
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Some  complicated  processes  occur  at  the  plasma-  metal  interface.  The  heat  release  at 
electrodes  under  Ohm  heating  of  their  surface  layer  and  under  absorption  of  radiation  from 
discharge,  as  well  as  a  result  of  electron  and  ion  plasma  heat  conduction  can  result  in  the 
electrode  melting  and  evaporation.  Nevertheless,  these  phenomena,  in  the  microsecond  range 
of  the  discharge  duration  even  in  high  power  APPT,  weakly  affect  APPT  parameters.  Indeed, 
it  has  experimentally  been  shown  that  the  losses  at  the  electrodes  do  not  exceed  10-15  %  in 
the  total  energy  balance.  The  boundary  layer  thickness  in  the  plasma  should  be  of  the  order  of 
a  few  particle  mean  free  paths  and  the  consideration  of  this  range  within  the  framework  of  die 
MHD-  model  does  not  seen  to  be  possible.  Therefore  it  is  expedient  to  formulate  the  boundary 
conditions  to  the  problem  to  exclude  the  boundary  layer  consideration  near  electrodes.  Within 
the  framework  of  the  MHD  -model  it  can  be  done,  assuming  that 


^(RH)=0 


Vr=o, 


At  the  electrode  surfaces  (when  R=Ri,  R=R2)  and 


(9) 


5vr  do  5H 
-^  =  0  ^  =  0  —  =  0 
dz  ^  dz  ^  dz  ^  (10) 

at  the  central  electrode  face  (when  z  =  zo,  R<Ri). 

As  for  the  velocity  Vz  at  z=zo,  R<Ri,  i.e.  at  the  cathode  cut,  two  versions  of  boundary 
conditions  are  possible: 

Vz  =  0,  the  cathode  end  is  impenetrable  for  the  plasma. 


C'Z  ,  (11) 

the  cathode  end  is  absolutely  penetrable  (divertor). 

At  the  APPT-axis,  when  R=0,  z>zo,  the  boundary  conditions,  follow  from  the  condition  of 
symmetry  for  all  the  values 


-^  =  0  —  =  0 

H=0,  VR=0,  ,  5R  .  (12) 

At  the  start  of  the  acceleration  channel  (z=0),  the  following  boundary  conditions  are  used 
(Ri=l): 


V,  =Vo,  v^  =  0,  p=Po,  T  =  To, 

R  .  (13) 

The  magnetic  field  at  the  point  (1,0)  is  determined  from  the  circuit  equation. 

Since  the  plasma  conductivity  is  a  finite  in  our  case,  one  should  formulate  the  boundary 
condition  for  the  magnetic  field.  At  the  output  from  the  APPT-  channel  (z=z0),  the  condition 
has  been  used. 
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It  is  often  used  in  the  MHD-  calculation.  It  corresponds  to  the  neglecting  of  a  radial  current 
density  component  at  the  acceleration  chaimel  output. 


III.2.6  Results  of  Simulation 


The  parameters  of  the  APPT  under  modeling  are  shown  in  Table  1. 
Table  1  Parameters  of  the  APPT  under  modeling 


Iimer  electrode  radius,Ri 

1  cm 

Radius  of  outer  electrode,  R2 

5  cm 

Length  of  irmer  electrode,  zo 

5-10  cm 

Length  of  outer  electrode,  z* 

17  cm 

Initial  voltage,  Vo 

1-3  kV 

Bank  capacitance,  Co 

300  pF 

Inductance,  Lo 

10’' H 

Average  mass  of  ion.  Mi 

16,6 

Temperature,  T 

1.5 -3  eV 

Input  plasma  density,  p 

10’"  -lOVcm" 

Input  plasma  velocity,  vo 

2-5  10"  cm/s 

Below  we  discuss  an  example  with  dimensionless  parameters  to  be:  Al=  50,  A2=3.6,  A3=l. 
Radial  distribution  of  ablated  mass  obtained  in  balance  experiments  is  in  inverse  proportion  to 
radius,  ~  1/r.  Total  ablated  mass  is  proportional  to  power  of  discharge.  For  relatively  high 
energy  flux  falling  onto  the  propellant  surface  propellant  flow  rate  is  proportional  to  this  flux. 
Approximation  of  inlet  condition  for  propellant  flow  rate  was  deducted  from  mass  loss 
measurements  and  modeling. 


ih  =  iho  +aR/R, — H^(R(„t) , 
dz 

riiQ  is  plasma  generated  by  ignition  device,  a-  is  the  coefficient  providing  a  correlation 
between  simulated  and  ablated  mass  loss,  Ri  -is  the  radius  of  central  electrode. 


As  it  was  already  marked,  the  discharge  in  APPT  is  accepted  conditionally  to  divide  into 
following  stages: 

a)  formation  of  the  discharge,  which  includes  a  temporary  interval  from  the  moment  of 
initiation  of  the  discharge  up  to  plasma  occurrence  at  the  outlet  of  the  accelerating  channel; 

b)  quasi-  steady  plasma  flow,  which  in  our  case  covers  an  interval  2t6.  In  this  stage  the 
distribution  of  parameters  of  plasma  in  the  accelerating  chaimel  varies  poorly; 

c)  damping,  when  a  current  and  power  of  the  discharge  are  small  in  comparison  with  their 
maximum  values. 

To  find  the  parameters  of  numerical  model  close  to  experimental  model,  firstly  the  current  - 
voltage  diagram  of  the  model  was  found  similar  to  experimental  one.  This  is  the  initial  point 
of  modeling.  Accordance  of  simulation  -  experimentation  in  integral  parameters  of  discharge 
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(current,-  voltage,  energy  transfer  to  APPT)  is  necessary  for  next  modeling  of  plasma 
parameters. 

The  current-voltage  diagram  of  APPT  discharge  under  consideration  is  shown  in  Figure  III.3. 


Figure  III.3.  Time  dependencies  of  a  current  and  hank  voltage  in  APPT  and  bank  voltage. 
to=  5  ns. 

The  maximal  current  and  forms  of  current  and  voltage  have  a  good  correlation  with 
experiment.  Figure  III.4,  Figure  III.5  give  the  distribution  of  plasma  density  and  magnetic 
field  in  the  acceleration  chaimel  for  moment  of  time  equal  1  (5ps ). 

These  distributions  are  changed  in  time  with  duration  as  current  and  voltage.  Other  words  at 
the  time  t  =  (1  -  2)  to,  quasi-steady  plasma  flow  takes  place.  Experiments  show  that  ionization 
layer  thickness  is  on  the  level  of  not  higher  than  0.1  mm.  So,  in  this  time  ionization  process 
does  not  effect  significantly  on  plasma  dynamics  near  the  inlet  of  the  thruster. 


H(R,Z) 


Outer 

electrode 


Figure  III.4.  Magnetic  field  distribution  in  APPT  accelerating  charmel  (t=  5  ps). 
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P(R.Z) 


Figure  III.5.  Plasma  density  distribution  in  APPT  accelerating  channel  (t=  5  |lis). 


Due  to  radial  distribution  of  plasma  density  in  the  in  the  inlet  to  the  acceleration  chaimel 
outlet  axial  velocity  of  plasma  have  weak  radial  dependence.  Other  picture  arises  for  radial 
velocity.  Inner  part  of  plasma  flow  is  focused  to  the  axis.  Plasma  located  near  outer  electrode 
has  negative  value  of  radial  velocity.  Accordingly  near  outer  electrode  (anode)  insignificant 
increase  in  plasma  density  arises.  This  does  not  change  noticeably  plasma  dynamics  in  the 
thruster. 

When  the  front  of  plasma  flow  reaches  the  end  of  the  central  electrode,  the  flow  with 
characteristic  tubular  structure  is  formed  in  this  area.  The  similar  structure  of  a  flow  exists  at 
the  end  of  the  central  electrode  and  further  the  distribution  of  plasma  parameters  varies 
poorly.  The  tubular  distribution  of  density  received  at  the  simulation,  has  a  good  correlation 
with  interferometric  results. 

Part  of  plasma  focused  in  the  axis  producing  plasma  focus  (the  area  of  increased  plasma 
density). 
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Figure  III.6.  Axial  velocity  distribution  in  APPT  accelerating  channel  (t=  5  (xs). 


Vr(R,Z)  t=1 


Figure  III.7.  Radial  velocity  distribution  in  APPT  accelerating  channel  (t=  5  |xs). 
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P(R.Z) 


Figure  IIL8.  Plasma  density  distribution  in  APPT  accelerating  channel  (t=  10  )is). 

Plasma  density  distribution  near  the  zero  of  a  current  (t=  10  jis)  shows  that  significant  part  of 
plasma  having  low  velocity  located  in  the  inlet  of  the  accelerating  channel.  Seems,  this  part  of 
a  plasma  at  the  end  of  first  half-period  of  a  current  create  late  time  ablation  mass  flow  visible 
in  the  interferogram. 

There  are  not  effectible  current  loops  in  volume  of  the  accelerating  channel  that  is  due  to 
rather  weak  dependence  of  plasma  velocity,  from  radius  and  large  size  of  average  plasma 
velocity.  Thus  the  skin-effect,  causing  occurrence  of  current  loops,  plays  a  smaller  role. 

As  a  whole  demonstration  of  two-dimensional  model  of  APPT  shows  that  this  adequately 
describes  formation  and  acceleration  of  a  plasma  flow  in  coaxial  accelerating  channel.  The 
simulated  parameters  of  considered  APPT  have  appeared  close  to,  measured  in  experiments. 
The  model  gives  a  correct  qualitative  picture  of  separate  stages  of  discharge  in  APPT: 

-plasma  flow  formation, 

-appearance  of  tubular  structure  at  the  end  of  the  central  electrode, 

-focusing  of  a  flow  etc. 

-appearance  of  low  accelerated  mass  near  the  zero  of  the  bank  voltage. 

The  features  of  propellant  flow  rate  determine  monotonous  character  of  distribution  of  density 
and  especially  axial  velocity  of  plasma,  and  magnetic  field  in  accelerating  chatmel. 

The  developed  numerical  model  of  APPT  is  a  tool  to  improve  thruster  parameters. 

III.3  Heat  Flux  Transfer  to  the  Propellant  and  Ablation  Losses. 

IIL3.1  Experimental  Stand  for  the  Propellant  Study 

The  quasi-steady  APPT  is  housed  in  a  rectangular  stainless  steel  tank  of  volume  near  1.2  m^. 
Prior  to  thruster  ignition,  the  tank  is  evacuated  down  to  approximately  3x10'^  Torr  by  an  oil 
diffusion  pump  and  mechanical  pump.  The  evacuation  velocity  is  equal  5000  1/s  under  the 
pressure  10^  Torr.  The  pulsed  power  supply  was  assembled  of  capacitors  with  the  total 
capacitance  up  to  300  pF.  Maximal  initial  voltage  is  up  to  5  kV.  The  capacitor  bank  is 
connected  with  the  feeder  through  the  low  inductive  cable  bridge  composed  of  48  high  current 
coaxial  cables.  The  estimated  cable  bridge  inductance  is  3  10’^  He.  The  inductance  of  vacuum 
feeder  is  5  10'*  H.  So,  bank  and  feeder  can  simulate  pulsed  power  source  of  APPT  in  wide 
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range  of  circuit  inductance  and  capacitance.  The  thruster  (Figure  III.9)  consists  of  a 
cylindrical  copper  anode  and  a  tungsten  cathode.  The  accelerating  coaxial  channel  is  near  9 
cm  long  and  10  cm  and  2  cm  in  diameter  respectively.  Outer  electrode  is  anode.  Teflon  is  the 
solid  propellant.  The  main  high  current  discharge  in  accelerating  coaxial  chaimel  is  triggering 
by  low  power  breakdown  plasma  moving  through  the  holes  in  the  inner  electrode  to  the  outer 
electrode 


a  b 

Figure  III.9.  Coaxial  APPT.  a-  face  view,  b  -  photo  of  the  discharge  in  the  APPT. 

in.3.1.1  Diagnostics 

Electrical  measurements  include  the  discharge  current,  voltage  and  the  magnetic  field 
measurements.  The  discharge  current  was  measured  by  the  Rogowski  coil  connected  to  low 
inductance  and  low  resistance  shunt.  The  voltage  across  the  capacitor  bank  and  at  other  points 
was  measured  by  capacitive  voltage  dividers.  The  high-speed  photography  of  the  discharge  is 
produced  with  a  high-speed  photo  camera.  Maximal  time  resolution  of  the  camera  in  streak 
mode  is  2  10‘*c. 

Energy  flux  propagated  to  propellant  surface  was  measured  with  the  tool  based  on  low 
dimensions  (1  mm  scale)  thermistor.  Sensor  is  located  into  hole  of  the  propellant  bar.  It  is 
glued  to  diin  mica  plate,  serving  as  a  support.  A  layout  of  involved  elements  is  shown  in 
Figure  III.  10.  Energy  flux  through  small  hole  in  the  propellant  bar  impacted  with  the  sensor. 
So,  sensor  stored  integral  flux,  which  produced  temperature  increase  and  a  change  in  the 
resistance.  A  comparison  bridge  circuit  is  used  to  measure  resistance  change  due  to 
temperature  increase.  Temperature  increase  occurs  after  near  1  s  after  shot.  Heat  transfer  from 
the  sensor  to  the  outer  space  has  typical  time  of  1  minute.  Sensor  is  adjusted  on  azimuthal 
direction  for  maximal  energy  flux  signal. 


Figure  III.  10.  Draft  of  the  measuring  tool  located  into  Teflon  insert  (a)  and  propellant 
bar  with  Teflon  insert  (b). 

Energy  flux  can  be  transferred  from  discharge  to  the  propellant  by  particles  and  by  radiation. 
To  outline  the  energy  flux  content  a  number  of  filters  were  used.  Data  on  filters  are  shown  in 
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Table  1.  Absorbance  dependencies  on  radiation  spectrum  are  shown  in  Figure  III.l  1.  Filters 
are  placed  on  the  surface  of  propellant  bar  closely  to  the  hole. 


Table  1.  Used  filters 


Material 

Thickness,  p,m 

Teflon 

4,  10,20 

Polyether 

6 

Nitrocellulose 

0.1 

Figure  III.l  1.  Absorbance  of  used  filters.  Absorbance  of  Teflon  and  polyether  is  given 
from  absorbance  of  Nitrocellulose  was  accounted. 


III.3.2  Results  of  Experiments 


APPT  discharge  under  consideration  has  not  symmetric  current  distribution.  (Figure  III.9-b). 
Such  a  discharge  is  similar  as  in  breach  fed  rail  geometry  thrusters.  Results  obtained  for  this 
discharge  include  current  and  voltage  distributions,  light  intensity  registration  fi:om  fiber 
inserted  near  energy  flux  registration  point.  Teflon  mass  loss  measurements  were  performed 
in  full  range  of  bank  energies.  Teflon  mass  loss  per  single  pulse  was  obtained  for  not  less  than 
300  pulses  and  each  stored  energy.  Time  dependency  of  the  discharge  parameters  are  shown 
in  Figure  III.  12. 
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Figure  III.  12.  Time  dependencies  of  the  discharge  current  I,  voltage  on  the  bank,  U  and 
radiation  flux  density,  H  (visible  wave  range,  optical  fiber  signal,  ab.  units). 


Calibration  of  energy  flux  measurements  was  produced  by  introducing  to  the  sensor  known 
quantity  of  heat.  It  is  realized  by  means  of  electrical  discharge  through  thin  platinum  foil 
glued  on  thermistor.  Calibration  signal  was  similar  working  one.  Typical  voltage  -  time 
dependency  of  the  sensor  for  1 80  J  bank  charge  is  shown  in  Figure  III.  1 3 .  Layout  of  sensors  is 
given  in  Figure  III.  14. 


Figure  III.  13.  Sensor  voltage  for  180  J  discharge 
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Figure  HI.  14.  Placement  of  sensors,  face  view 


Radiation  flux  determined  with  polyether  and  nitrocellulose  filter  was  not  more  than  1%  of 
energy  flux  measured  without  filters.  Installation  of  6  fxm,  10  pm,  and  20  pm  Teflon  filters 
produced  the  signal  at  level  5%  of  total  (particle  +radiation).  Teflon  filters  were  not 
destroyed.  Nitrocellulose  filters  were  usually  destroyed  in  the  discharge  after  one  shot. 
Measured  energy  fluxes  are  shown  in  Figure  111.15.  Curve  “1”  show  energy  flux 
measurements  from  centered  sensor  position  (Figure  III.  14,  #1).  In  this  case  all  particle  and 
radiation  flux  should  impact  with  sensor.  Curves  “2”and  “3”  obtained  with  radial 
displacement  of  the  sensor.  “2”  -  surface  of  the  sensor  located  in  1  mm  from  side  of  hole.  “3” 
-  surface  of  the  sensor  located  in  3  mm  from  side  of  hole.  In  position  “3”  sensor  does  not  see 
area  of  plasma.  So,  as  seen  from  Figure  111.15,  conductive  energy  transfer  to  the  propellant 
exceed  10%  of  total  transferred  energy.  There  are  no  visible  damages  on  the  sensor  surface 
during  measurements. 


Figure  III.  15.  Measured  energy  fluxes  propagated  onto  propellant  (1-4)  and  mass  loss  in 
dependence  on  the  bank  energy.  1-  without  filter  (centered);  2-  without  filter,  2  mm 
displacement  from  hole,  3  -  without  filter,  3  mm  displacement  from  hole;  4-6  pm  Teflon 
film  filter. 

The  electron  density  of  plasma  in  the  near  propellant  surface  area  was  measured  with  Mach- 
Zehnder  interferometer.  The  ruby  laser  -  in  a  quasi-continuous  generation  mode  with 
wavelength  0.69  p  -  was  used  as  an  illuminator.  The  propellant  plane  was  aligned  in  parallel 
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to  the  optical  interferometer  axis  with  the  accuracy  up  to  3’.  The  interference  pattern  shift  was 
registered  with  a  streak  camera.  Streak  photo  interferogram  of  near  propellant  area  are  given 
in  Figure  III.16.The  time  resolution  was  5  10'*s.  Electron  density  depended  on  the  input 
voltage  and  measured  at  the  level  of  10*^cm"^  for  0.5  kV  (inter-electrode  gap  voltage)  with 
some  deviation  due  to  shot-to-shot  changes  in  thruster  firing. 


Figure  III.  16.  Streak  photo  interferogram  of  near  propellant  area.  Time  of  streak  is  1.5  ps. 
Estimated  electron  density  is  lO'^cm'^.  Voltage  on  electrodes  is  near  0.5  kV  (4  cm  inter¬ 
electrode  gap). 

It  is  seen  from  Figure  III.  16,  that  there  are  oscillations  in  electron  density  at  a  distance  -  near 
(0.2  -0.3)  mm  from  propellant  surface.  Estimated  period  of  these  oscillations  (instabilities)  is 
70  ns. 

Dependence  of  mass  loss  distribution  in  the  inlet  of  the  accelerating  channel  is  necessary  for 
semi-phenomenological  modeling.  Model  needs  the  azimuthally  symmetric  conditions.  To 
initiate  azimuthally  symmetric  discharge  a  small  quantity  of  Aluminum  (significantly  smaller 
than  propellant  loss)  was  evaporated  and  deposited  onto  propellant  surface.  Accordingly 
vacuum  switch  was  used  in  the  electric  circuit.  The  discharge  became  azimuthally  symmetric. 
An  example  of  radial  distribution  of  a  mass  loss  in  APPT  obtained  from  balance 
measurements  is  shown  in  Figure  III.17. 


Figure  III.  17.  Dependence  of  mass  loss  distribution  in  APPT  in  azimuthally-symmetric 
discharge.  Bank  energy  is  180  J. 
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III.3.2.1  Measurements  of  Energy  Dissipated  in  a  Thermal  Skin-Layer. 

To  measure  energy  dissipated  in  a  thermal  skin-layer  the  discharge  was  produced  in  the 
surface  of  thin  Teflon  film,  which  was  closely  pushed  to  the  sensor  as  shown  in  Figure  III.  18. 


Energy  flux 


iiiiiiiiniiiiiiiiiuiiuiiiiiui 


Figure  III.  1 8.  Method  of  thermal  energy  measurement  in  the  skin-layer 


Measured  energy  bit  dissipated  in  Teflon  propellant  per  one  discharge  is  shown  in  Figure 

ni.l9. 


Bank  energy,  J 


Figure  III.  19.  Thermal  energy  bit  dissipated  in  Teflon  propellant  at  the  discharge. 


III.3.2.2  Measurements  of  Radiative  Part  in  the  Energy  Flux . 

Measurements  of  radiative  flux  with  the  use  of  available  filters  are  limited  of  wavelength  near 
160  nm.  It  is  shown  above  that  radiative  part  of  energy  flux  with  wavelength  more  than  160 
nm  does  not  exceed  0.05  of  total  measured  energy  flux.  Shorter  radiation  can  be  studied  using 
mirrors.  Tantalum  mirrors  adjusted  to  sensor  were  used  to  detect  radiative  part  of  energy  flux. 
Reflectivity  of  Tantalum  in  near  100  nm  for  normal  slope  is  near  20%. 

Adjustment  and  de-adjustment  of  reflecting  mirror  was  produced  for  transparent  source 
illumination  located  near  propellant  surface.  With  adjusted  mirror  sensor  should  detect  optical 
signal.  With  de-adjusted  mirror  optical  signal  did  not  fall  onto  sensor  (Figure  III.20,  a,  b) 
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Figure  III.20.  Determining  of  energy  flux  structure;  a-  mirror  adjusted  to  sensor  for  measuring 
of  radiation,  b-mirror  does  not  register  radiation  (not  adjusted). 


Taking  into  account  some  uncertainties  due  to  shot-to-shot  variations  in  thruster  firing,  sensor 
did  not  detect  appreciable  energy  flux  transferred  by  radiation  (Figure  III.21).  So,  noticeable 
addition  associated  with  radiative  energy  flux  has  not  detected. 


0  100  200  300  400  500  600 

Bank  energy,  J 


Figure  III.21.  Determining  of  energy  flux  structure,  results  of  measuring.  Signals  of 
sensor  for  adjusted  and  de-adjusted  mirrors. 
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In  the  APPT,  the  considerable  energy  fraction  released  on  the  propellant  bar  surface  may  be 
transferred  by  particles  and  ultraviolet  radiation  in  the  wavelength  range  100  nm.  The 
absorption  depth  of  this  radiation  in  polymers  does  not  exceed  10'^  cm  meanwhile  the 
thickness  of  an  propellant  material  heating  due  to  the  heat  conduction  is  about  10^  cm. 
Schematic  of  the  energy  transfer  in  APPT  is  shown  in  Figure  III.22.  At  present,  die  relative 
role  of  conductive  and  radiation  mechanisms  for  energy  transfer  to  the  propellant  is  not 
exactly  determined.  Now  conductive  mechanism  is  more  expanded  and  developed 
Probably  the  main  mechanism  of  energy  transfer  to  propellant  depends  on  propellant  and 
thruster  operation  mode. 

In  considered  experiment  we  measure  energy  flux  through  a  low  dimensions  hole  in 
propellant  bar.  Most  probably  this  hole  destroys  boundary  layer  near  the  propellant  surface, 
because  its  thickness  exceeds  significantly  the  thickness  of  boundaiy  layer.  In  this  case  we 
possibly  measure  energy  flux  falls  on  boundary  layer.  For  this  reason  experiments  with  use  of 
Teflon  filters  are  more  reliable. 

Boundary  layer 


Anode 


Figure  III.22.  Schematic  of  the  energy  transfer  in  APPT 


The  approximation  of  energy  flux  transferred  onto  the  propellant  is  shown  in  Figure  III.23. 
This  flux  could  be  used  as  a  boundaiy  eondition  for  accounting.  Accounts  of  Teflon 
degradation  based  on  thermal  degradation  kinetics  (Annex  1)  have  been  done.  The  results  of 
accounts  are  shown  in  the  Table  2:  Hmax  -maximal  power  of  energy  flux,  q-  total  energy  flux 
per  cm^,  ablated  mass  Maw,  Qskin-  energy  of  thermal  skin  layer,  Qabi,-  Teflon  destruction 
energy  and  Ts, -maximal  temperature  of  the  propellant  surface. 
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Figure  111.23.  Time  dependence  of  energy  flux  onto  propellant. 


Table  2.  Results  of  accountings 


Hmax, 

kW/cm^ 

q  =  jHdt 
J/cm^ 

Mabl, 

pg/cm 

Qskiii) 

J/cm^ 

Qabl» 

J/cm^ 

Ts,®K 

12.2 

0.11 

0 

0.126 

0 

705 

24.3 

0.218 

0.75 

0.224 

3.10'^ 

1090 

48.6 

0.435 

39 

0.260 

0.168 

1170 

73.2 

0.653 

83.3 

0.262 

0.375 

1190 

97.6 

0.871 

130 

0.258 

0.603 

1210 

122 

1.09 

184 

0.252 

0.862 

1230 

The  results  of  accounting  show  that  Teflon  begins  to  destruct  from  surface  temperature  > 
llOO^K.  Further,  temperature  increases  insignificantly  with  energy  flux  increasing.  When  the 
power  of  energy  flux  increases  8  times,  tiie  temperature  increases  on  100®.  For  Ts=  1200®K 
and  for  energy  flux,  equal  shown  in  Figure  III.23,  thermal  skin-layer  contain  0.260  J.  High 
mass  loss  case  is  demonstrated  in  Figure  III.24,  where  H,  Mabl,  Ts  and  e  are  shown  for 
Hmax=73  kW/cm^.  Temperature  distribution  and  Teflon  destruction  rate  for 
Hmax=73  kW/cm^  (t=4  ps)  are  shown  in  Figure  III.25. 
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Figure  111.25.  Temperature  distribution  and  Teflon  destruction  rate  for  Hmax=73  kW/cm' 
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Figure  III.26.Dependencies  of  Hmax  h  Qskin  of  ablated  propellant  mass. 


Maximal  value  of  energy  flux  power  to  the  propellant  Hmax.  necessary  to  evaporate  mass 
Mabl,  In  dependence  of  Mabl  is  shown  in  Figure  in.26  together  with  thermal  energy  in  skin- 
layer.  It  follows  from  the  last  figure  that  Teflon  degradation  for  measured  energy  fluxes 
falling  onto  the  propellant  can  take  place  if  temporal  or  spatial  inhomogeneities  in  energy  flux 
occur.  In  this  case  measured  energy  is  average  value  and  can  be  significantly  lower  than  local 
values  of  energy  flux. 

Such  model  can  explain  the  mass  loss  in  a  frame  of  thermal  degradation  kinetics  if  maximum 
energy  flux  to  the  propellant  exceeds  25  kW/cm^  (Integral  flux  more  than  0.2  J/cm^).  So,  in  a 
frame  of  used  thermal  degradation  model  of  Teflon,  for  measured  energies  shown  in  Figure 
III.  15  and  Figure  HI.  19  degradation  of  Teflon  is  doubtful.  However  inhomogeneities  in 
energy  flux  distribution  can  explain  observed  mass  loss  for  relatively  low  energy  fluxes. 
Figure  III.  16,  where  oscillations  in  electron  density  at  a  distance  -  near  (0.2  -0.3)  mm  from 
propellant  surface  are  shown  could  confirm  inhomogeneities  in  energy  flux  distribution. 
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II1.4  Study  of  plasma  parameters  and  discharge  optimization 

111.4.1  Study  of  the  charged  plasma  component  concentration  in 

discharge  channel 


The  electron  density  is  the  fundamental  parameter  of  APPT  plasma.  The  distinctive  features 
of  APPT  physical  processes,  current  distribution,  acceleration  mode  of  plasma  blob  and  thrust 
performances  such  as  impulse  bit,  plasma  velocity  (specific  thrust),  thrust  efficiency  etc  are  in 
many  respects  dependent  on  the  value  and  the  distribution  of  electron  density  as  well  as  ion 
density. 

A  method  of  plasma  laser  interferometry  was  used  for  measuring  the  electron  density  inside 
the  discharge  channel  of  next  generation  APPT  and  in  the  plasma  blob  outside  die  discharge 
channel.  The  advantage  of  this  method  lies  in  the  fact  that  plasma  under  study  is  not  disturbed 
and  that  this  method  allows  spatial  and  temporal  distributions  of  electron  density  be  measured 
in  non-equilibrium  plasma  with  acceptable  accuracy. 

Test  data  on  the  measured  electron  density  in  the  discharge  chaimel  zone  between  the  Teflon 
bars  were  obtained  while  fulfilling  this  work  for  the  first  time.  Plasma  flow  is  formed  in  this 
zone  as  a  result  of  initial  breakdown  of  the  discharge  gap,  propellant  ablation  and  ionization, 
and  a  current  sheet  is  formed  also  that  is  later  accelerated  by  electromagnetic  forces.  Thus, 
studies  for  the  electron  density  distribution  in  the  zone  between  the  bars  are  very  important  for 
the  understanding  of  operating  processes  in  APPT. 


III.4.2  Measurement  procedure 

Laboratory  model  of  next  generation  APPT  with  the  stored  energy  of  about  100  J  was  used  in 
the  tests  on  measuring  the  electron  density  distribution  and  its  variation  in  time.  Main 
elements  of  the  discharge  channel  unit  of  that  model  are  shown  in  Figure  III.27.  Discharge 
chatmel  comprises  electrodes  1  and  2,  two  propellant  (Teflon)  bars  3  and  the  end  insulator  4. 
A  side  location  of  Teflon  bars  was  used  in  the  APPT  model.  The  anode  is  installed  in  die  base 
5  and  cathode  -  in  the  cover  8.  The  APPT  electrodes  are  manufactured  of  copper.  The  end 
insulator  is  made  of  alundum  plates,  which  have  been  pasted  together.  They  are  mounted  in 
yoke  6  made  of  glass-cloth  laminate.  The  Teflon  bars  3  are  installed  in  yoke  6  closely  to  the 
end  insulator  4.  Distance  betwfeen  electrodes  of  70  mm  in  length  was  38  mm,  while  the 
longitudinal  size  of  Teflon  bars  was  30  mm. 

The  APPT  discharge  initiation  is  carried  out  by  three-electrode  igniter  7  installed  in  cover  8. 
The  igniter  is  powered  from  the  discharge  initiation  unit.  Capacitor  bank  consisted  of  four 
low-inductive  capacitors. 

APPT  model  was  tested  in  the  vacuum  tank  at  the  residual  gas  pressure  of  about  10^  torr.  The 
known  procedures  and  measuring  equipment  described  in  were  used  for  measuring  the 
APPT  main  parameters  (current  pulse,  energy,  thrust,  specific  thrust  impulse,  mean  propellant 
flow  rate  per  pulse,  thrust  efficiency,  etc.). 

Electron  density  Ne  in  the  APPT  plasma  was  measured  by  laser  interferometer  with 
photoelectric  signal  registration  In  preparation  for  experiments  the  interferometer  was 
updated  with  the  aim  to  improve  the  spatial  and  temporal  resolution  and  to  increase  the 
measurement  accuracy.  The  laser  power  of  updated  interferometer  was  10  times  greater  in 
comparison  with  the  laser  used  before  The  noise-immunity  for  electrical  circuit  was 
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improved  and  the  damping  system  was  refined  for  reducing  the  mechanical  vibrations.  The 
up^ted  laser  interferometer  provides  for  measuring  the  electron  density  in  the  range  from 

3T0'^  to  lO'®  cm'^  with  the  spatial  resolution  of  ~  2  mm  and  temporal  resolution  of  ~  1  ps. 

The  optical  scheme  of  interferometric  measurements  is  presented  in  Figure  III.28.  The 
updated  interferometer  used  during  the  experiments  was  one  of  the  versions  of  classical  three- 
mirror  laser  interferometer, Active  resonator  of  the  interferometer  is  file  15  mW  single¬ 
mode  helium-neon  laser  1.  Passive  resonator,  in  which  the  investigated  APPT  plasma  is 
placed,  is  formed  by  an  output  mirror  of  the  laser  and  external  spherical  mirror  5.  The  mirror 
5  of  this  resonator  with  the  focal  distance  of  1  m  has  the  reflection  coefficient  of  95  %  for  the 
wavelength  X  =  6328  A. 

To  generate  periodic  oscillations  of  optical  length  in  the  passive  resonator,  the  external  mirror 
is  fixed  to  the  piezoelectric  ceramic  cylinder  providing  periodic  shift  for  the  mirror  with  a 
frequency  of  50  Hz  when  the  modulating  voltage  of  »10  V  is  applied.  The  movable 
diaphragm  6  is  used  while  choosing  an  interference  order,  which  is  necessary  for 
measurements.  The  input  and  output  of  laser  radiation  through  the  vacuum  chamber  is  carried 
out  by  using  quartz  windows  of  5  cm  in  diameter. 

The  registration  of  optical  signal  is  carried  out  by  photomultiplier  8  and  digital  memory 
oscilloscope.  The  highest  spectral  sensitivity  of  photomultiplier  corresponds  to  the  wave 
length  of  ~  600  nm  that  is  rather  close  to  the  laser  wave  length  (X  =  632.8  ran). 

During  the  experiments  the  optical  scheme  was  adjusted  so  that  the  amplitude  of  signal 
measured  by  the  photomultiplier  did  not  exceed  500  mV,  that  corresponded  to  linearity 
boundary  of  this  device.  The  number  of  laser  ray  passes  in  the  passive  resonator  is  evaluated 
by  the  value  of  10. ..  12. 

The  photomultiplier  was  powered  by  a  high-voltage  stabilized  power  supply.  The  nominal 
value  of  voltage  was  equal  to  1050  V.  A  resolution  of  the  photomultiplier  was  estimated  by 
the  value  of  ~  1  ps.  This  may  be  apparently  considered  as  satisfactory  because,  on  the  one 
hand,  this  value  essentially  exceeds  relaxation  time  of  physical  processes  in  resonators  of  die 
interferometer  and,  on  the  other  hand,  it  is  much  less  than  the  duration  of  discharge  process  in 
APPT(~  15...20ps). 

The  noise-immunity  of  the  electrical  circuit  was  ensured  by  grounding  the  experimental 
instrumentation.  Besides,  all  high-voltage  and  a  part  of  low-voltage  circuits  were  passed 
through  ferrite  rings,  that  allowed  the  reduction  of  noise  down  to  an  acceptable  level. 

Some  precautionary  measures  against  light  interferences  in  the  receiving  part  of 
photomultiplier  accompanying  the  APPT  operation  were  taken  also.  At  the  exit  of  the  vacuum 
chamber  the  light  screen  with  a  hole  of  10  mm  for  a  ray  passing,  neutral  light  filters  and 
movable  diaphragm  of  0.7  mm  in  diameter  were  installed.  Besides,  a  photomultiplier  was 
located  at  distance  of  ~  Im  from  APPT  model,  and  it  was  also  housed  in  a  package,  which 
was  opaque  to  light  and  had  an  inlet  opening  of  ~  2  mm  for  laser  radiation  passing.  All 
measurements  were  conducted  at  exterior  black-out. 
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Figure  111.27.  Discharge  channel  unit  schematic 

1  -  cathode,  2  -  anode,  3  -  teflon  bars,  4  -  back  insulator, 
5  -  basis,  6  -  yoke,  7  -  igniter,  8  -  cover. 


Figure  III.28.  Optical  scheme  for  interferometric  measurements. 

1  -  He-Ne  laser,  2  -  diaphragm,  3  -  vacuum  chamber,  4  -  APPT  model,  5  -  mirror,  6  - 
diaphragm,  7  -  filter,8  -  photomultiplier. 


To  reduce  a  level  of  optical  system  vibrations,  all  devices  of  the  laser  interferometer  (helium- 
neon  laser,  mirrors,  diaphragms,  filters,  photomultiplier)  are  located  on  a  metal  plate  with  the 
mass  of  more  than  100  kg  suspended  by  four  rubber  braids  from  the  metal  fi'ame.  This  frame 
surrounding  the  vacuum  chamber  is  mounted  rigidly  to  a  metal  construction  of  200  kg  in  mass 
placed  on  the  floor  of  experimental  stand  by  using  the  multi-layer  damping  spacers  made  of 
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elastic  materials.  After  each  cycle  of  measurements  the  plate  with  laser  interferometer  was 
arrested  by  a  special  device,  that  allowed  the  off-loading  of  rubber  suspending  braids.  In 
general,  the  whole  constmction  has  ensured  the  reduction  in  level  of  optical  system  vibrations 
by  the  factor  of  1 5 ..  .20. 

Thus,  during  experiments  protection  was  provided  for  the  electric  circuit  of  interferometer  and 
measuring  equipment  against  different  noises.  The  amplitude  of  noises  in  the  measuring 
circuit  (mechanical,  electromagnetic,  thermal)  comprised  less  than  10%  of  the  main  signal 
amplitude. 

The  synchronization  of  measurements  was  carried  out  by  multi-channel  pulse  generator 
making  it  possible  to  adjust  measurement  process  in  time  with  an  accuracy  of  0.1  p.s.  All 
incoming  data  was  recorded  by  two  digital  memory  oscilloscopes.  The  current  pulse  of  APPT 
model  and  the  interferogram  of  APPT  operating  process  were  recorded  by  one  oscilloscope 
and  the  instrumental  function  of  the  interferometer  was  fixed  by  the  other  oscilloscope  at  time 
of  APPT  model  switching.  Experimentally,  from  6  to  10  measurements  of  electron  density  Ne 
were  carried  out  in  each  given  spatial  point  and  then  measurement  results  were  averaged. 
Electron  density  measurement  in  plasma  by  application  of  laser  interferometer  is  based  on  the 
registration  of  changing  the  phase  of  laser  radiation  passing  through  plasma  under  study  [4]. 
This  change  of  the  phase  A(p  is  proportional  to  electron  density  Ne  and  plasma  size  d  in  the 
direction  of  laser  ray 

A^  =  — -lN^d  =  const -N^d,  (1) 

me 

where  X  is  the  wave  length  laser  radiation,  e  and  m  are,  correspondingly,  charge  and  mass  of 
electron,  c  is  the  light  velocity.  Thus,  the  information  on  the  electron  density  value  averaged 
in  the  direction  of  laser  ray  may  be  obtained  by  measuring  the  phase  difference  A(p  if  the 
plasma  size  d  is  known. 

A  method  of  laser  interferometiy  may  be  used  effectively  for  plasma  diagnostics  when  the 
following  requirements  are  met  [ 

1 .  Laser  radiation  should  not  be  absorbed  in  plasma. 

2.  Radiation  entering  plasma  should  keep  its  initial  direction,  and  laser  ray  deflection  at  the 
electron  density  gradients  should  be  negligibly  low. 

3.  Modulation  of  radiation  intensity  as  it  leaves  plasma  should  be  mainly  substantiated  by 
the  interference  effect,  and  the  effect  of  the  polarization  plane  rotation  for  the  laser 
radiation  (Faraday  effect)  upon  its  intensity  at  a  presence  of  magnetic  fields  in  plasma 
should  not  be  considerable. 

It  was  revealed  as  a  result  of  analytical  assessments  and  experimental  check  for  the  Ne 
gradient  presence  in  plasma  by  a  shade  method  that  for  the  plasma  with  parameters  equal  to 
the  APPT  under  study  the  above  requirements  are  met.  Besides,  assessment  by  formulae 
presented  in  showed  that  under  the  conditions  of  these  tests  the  plasma  refraction  is 
defined  by  electrons,  and  the  contribution  of  atoms  may  be  neglected.  Thus,  application  of 
laser  interferometry  for  measuring  the  electron  density  was  quite  correct  in  that  work. 

At  the  presence  of  APPT  plasma  in  the  passive  resonator  its  optical  length  varies  and 
according  to  this  the  fundamental  frequency  of  radiation  leaving  the  resonator  varies.  This 
frequency  shift  Av  is  defined  by  the  formula 

Av  =  (An-d)^,  (2) 

where  An  d  is  the  change  in  optical  length  caused  by  plasma. 
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An  is  the  change  of  refraction  index  in  the  resonator, 
d  is  the  size  of  plasma  along  the  laser  ray  direction, 

Lr  is  the  resonator  length, 

X  is  the  wave  length  of  laser  radiation. 

As  the  three-mirror  interferometer  operates  as  an  amplitude-frequency  discriminator  with 
linear  performance,  the  frequency  shift  caused  by  plasma  will  cause  proportional  change  in 
laser  power  AE  at  the  exit  of  the  interferometer 

Av 

AE  =  E,„ax—  ,  (3) 

o 


where  E^ax  is  the  maximum  laser  power, 

5  is  the  interval  of  frequency  variation. 


The  electron  density  in  plasma  is  determined  on  the  basis  of  interferometer  signals  recorded  at 
the  presence  of  plasma  in  the  passive  resonator.  Periodic  displacement  of  external  mirror  in 
the  passive  resonator  provides  the  registration  of  interferometer  instrumental  function  used  as 
a  calibration  signal.  The  typical  instrumental  function  obtained  during  one  of  the  experiments 
is  presented  in  Figure  III.29. 

Measurement  of  the  photomultiplier  signal  amplitude,  which  is  proportional  to  AE,  and 
measurement  of  instrumental  function  parameters  (amplitude  Gmax»  front  duration  tin  and 
period  Tm)  makes  possible  to  define  Q  value  as  a  product  of  electron  density  Ne  by  the  plasma 
size  d  in  the  direction  of  laser  ray  (Q  =  Ne-d)  characterizing  the  radiation  phase  change  at  the 
plasma  appearance  in  resonator.  Equation  for  Q  obtained  on  the  basis  of  known  relations  for 
the  plasma  refraction  factor  and  laser  power  change  has  the  following  form: 


T,.. 


where:  Emax  -  maximum  power  of  the  laser; 

AF  -  the  value  of  signal  under  study  received  by  photomultiplier; 

Gmax  -  signal  amplitude  for  the  instrumental  function; 

Atin  -  duration  of  the  instrumental  function  front  (Figure  III.29); 

Tin  -  instrumental  function  period  (Figure  in.29); 

X  -  longitudinal  coordinate  for  the  point  under  study  in  the  direction  of 
APPT  axis; 

y  -  transversal  coordinate  for  the  point  under  study  in  the  direction 
between  the  electrodes; 
t  -  time. 

Electron  density  Ne  was  measured  in  the  discharge  channel  in  the  zone  between  the  Teflon 
bars  and  in  the  interelectrode  gap  after  the  bar  cut,  as  well  as  outside  the  discharge  channel 
after  the  electrode  cut. 

A  diagram  for  the  location  of  lines  along  which  plasma  was  examined  by  ftie  laser  ray  while 
measuring  Ne  is  shown  in  Figure  III.30.  The  ray  of  helium-neon  laser  crosses  the  APPT 
plasma  under  study  in  parallel  to  the  electrode  plane,  i.e.  along  the  line  being  parallel  to  z  axis. 
Points  shown  in  Figure  III.30  are  the  points  in  which  lines  of  plasma  examination  intersect  the 
plane  being  perpendicular  to  the  electrodes  and  including  the  geometric  axis  of  APPT 
channel.  Ne  was  measured  on  the  discharge  chatmel  axis  (line  I  - 1),  near  the  caftiode  (line  11  - 
n),  and  near  the  anode  (line  III  -  III),  as  well  as  on  the  line  IV  -  IV  being  perpendicular  to  the 
direction  of  plasma  blob  motion  outside  the  APPT  at  a  distance  of  5  cm  from  the  electrode 
cut. 
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•  -  an  investigated  point 
P  -  the  direction  of  plasma  blob  motion 
1  -Teflon  bars 


Figiure  111.30.  Arrangement  of  points  for  interferometric  measurements  in  APPT 


One  hole  of  about  4  mm  in  diameter  was  made  in  each  Teflon  bar  in  accordance  with  the 
location  of  line  of  plasma  examination  for  the  laser  ray  to  pass  through  them.  These  holes 
were  consequently  closed  by  Teflon  caps  after  measurements  were  finished  in  each  of  the 
examination  directions.  According  to  the  data  of  comparative  tests,  a  presence  of  one  pair  of 
holes  in  the  propellant  bars  does  not  influence  practically  the  APPT  characteristics  (integral 
parameters  of  the  model  varied  by  ~1 . .  .2%). 

With  the  known  plasma  blob  extension  d  along  z  axis,  an  electron  density  averaged  along  the 
laser  ray  direction  may  be  defined  by  (4)  for  the  studied  points  of  APPT. 

Within  the  volume  between  Teflon  bars,  the  blob  size  d  along  z  axis  is  a  fixed  value  defined 
by  the  discharge  channel  width.  For  the  other  points  under  study  (Figure  ni.30),  d  was 
defined  similarly  to  the  procedure  presented  in  i.e.  by  measuring  the  blob  boundaries 
using  the  traces  left  by  plasma  blob  on  the  dielectric  screen  during  the  APPT  model  operation 
in  the  frequency  mode.  This  screen  was  mounted  perpendicularly  to  the  direction  of  blob 
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motion  at  different  distances  from  the  electrode  cut.  The  variation  of  plasma  blob  extension 
d(x,0)  along  x  axis  in  the  plane  being  parallel  to  electrodes  and  including  the  channel  axis 
while  plasma  blob  moves  after  the  electrode  cut  is  presented  in  Figure  III.3 1 . 

In  addition  to  the  above  method  and  for  the  verification  of  test  results  with  screens,  plasma 
size  d  was  also  determined  by  interferometric  measurements.  The  plasma  extension  d  was 
assessed  by  the  interferometer  signal  decay  and  thus  by  the  reduction  of  Q(z)  parameter  when 
plasma  was  examined  in  the  direction  parallel  to  the  plane  of  electrodes  (y  axis). 

During  that  test,  the  APPT  model  was  turned  by  90*^  relative  to  the  longitudinal  axis  x  for  the 
laser  ray  to  be  directed  in  parallel  to  y  axis  (Figure  III.30).  Distribution  of  the  parameter  Q  = 
Ne-d  characterizing  the  variation  of  phase  for  the  laser  radiation  was  defined  by  interferometer 
along  z  axis.  Measurements  were  made  with  5  mm  step  in  the  cross-section  of  discharge 
channel  located  in  the  middle  of  the  electrode  part  after  the  Teflon  bar  cut,  and  in  the  cross- 
section  of  the  channel  at  the  electrode  cut.  The  obtained  Q(z)  distributions  allow  the  extension 
of  plasma  blob  d(z)  be  assessed  by  defining  the  blob  boundaries  on  the  basis  of  Q  value  drop 
down  to  the  given  level  with  the  growth  of  distance  from  the  APPT  channel  axis. 
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Figure  111.32  illustrates  the  variation  of  value  Q  =  Ne-d  in  the  areas  adjacent  to  the  electrodes, 
in  the  channel  cross-section  located  in  the  middle  of  the  electrode  part  after  the  Teflon  bar  cut. 
Points,  in  which  Q  value  dropped  five  times  comparing  to  the  maximum  value,  were 
considered  as  the  boundaries  of  plasma  blob  along  the  examination  line  while  the  test  data 
was  processed.  Blob  size  d  determined  by  measuring  the  distributions  Q(zj  differed  from  d 
size  determined  as  a  result  of  tests  with  screens  by  no  more  than  20%.  Such  agreement  of  test 
results  is  quite  satisfactory. 


213 


ISTC 


FINAL  TECHNICAL  REPORT 


Partner  Project  U231 7p 


d,  mm 


! - ! 

electrode 


Figure  IIL31.  Variation  of  plasma  blob  size  along  the  APPT  axis 
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Figure  IIL32  Variation  of  Q=Ne  d  parameter  in  areas  adjacent  to  electrodes 

IIL4.3  Test  results  and  their  discussion 

Integral  characteristics  of  next  generation  APPT  model  were  studied  within  the  stored 
energy  range  of  70. ..  100  J  at  the  first  stage  of  tests.  Main  parameters  of  this  model,  averaged 
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over  the  results  of  five  measurements,  are  presented  in  Table  S.Parameters  of  the  APPT  model 
under  study  are  substantially  better  than  parameters  of  the  typical  APPT’s  .  For  comparison  it 
should  be  pointed  out  that  in  typical  APPT  at  energy  W  =  100  J  the  thrust  efficiency  T|t 
amounted  to  0. 12,  impulse  bit  Pbh  =  2.2  mN-s  and  specific  impulse  Isp  =  1 1  km/s  [2], 


Table  3  Characteristics  of  the  next  generation  APPT  model 


Discharge  energy,  J 

70 

100 

Thrust  impulse  bit,  mN  s 

1.82 

2.7 

Teflon  consumption  per  pulse,  mg 

0.12 

0.15 

Thrust  efficiency 

0.20 

0.24 

Specific  impulse,  km/s 

15.4 

18 

The  discharge  current  oscillogram  for  the  APPT  model  under  study  at  the  energy  of  100  J 
is  presented  in  Figure  III.33.  It  is  seen  that  the  discharge  in  the  chaimel  of  next  generation 
APPT  is  close  to  aperiodic.  Current  amplitude  was  27.9  kA  during  the  first  half-period  of  the 
discharge  and  5.4  kA  during  the  second  half-period  (at  the  opposite  current  direction).  The 
duration  of  the  discharge  first  half-period  was  8.9  ps  that  was  close  to  the  duration  of  the 
process  of  plasma  blob  formation.  TTiis  provided  optimum  conditions  for  the  electromagnetic 
acceleration  of  plasma  ^^'^Electron  density  Ne  in  the  discharge  channel  between  the  Teflon 
bars  was  measured  for  the  first  time  while  studying  the  parameters  of  next  generation  APPT 
by  the  method  of  laser  interferometry.  Ne  measurements  in  the  points  located  in  the 
interelectrode  gap  after  the  bar  cut  and  outside  the  discharge  charaiel  (Figure  III.30),  which 
had  been  presented  earlier  in  were  made  again  in  order  to  obtain  data  on  the  Ne  variation 
along  the  entire  length  of  discharge  channel  during  one  run  of  measurements.  It  follows  from 
the  analysis  of  data  of  interferometric  measurements  obtained  inside  the  discharge  channel 
that  the  duration  of  interferometric  signal  corresponds  to  the  duration  of  half-period  for  the 
discharge  current.  This  may  probably  witness  of  the  fact  that  ablation  processes  inside  the 
channel  proceed  during  the  time  interval  of  7-8  ps  from  the  discharge  start  mainly. 
Interferometric  measurements  of  the  electron  density  Ne  showed  that  the  current  sheet  being 
formed  inside  the  discharge  chaimel  has  a  longitudinal  size  (along  x  axis)  of  5-7  cm,  and  that 
there  is  a  plasma  area  in  the  central  part  of  the  current  sheet,  in  which  the  electron  density  is 
as  high  as  the  maximum  value. 

Time-dependent  distributions  of  maximum  electron  density  Ne(x,  0,  t)  along  the  next 
generation  APPT  axis  while  plasma  moves  in  the  zone  between  Teflon  bars,  in  the 
interelectrode  gap  after  the  bar  cut  and  outside  the  discharge  channel  is  shown  in  Figure 
III.34.  Time  moments  since  the  discharge  start,  which  correspond  to  the  appearance  of  plasma 
layer  with  maximum  values  of  electron  density  in  the  APPT  points  under  study,  are  presented 
near  the  measured  values  of  Ne  also. 

The  trend  of  Ne  variation  during  the  plasma  motion  through  the  space  between  bars 
corresponds  to  the  variation  of  discharge  current  J(t)  for  the  next  generation  APPT,  in  which 
the  duration  of  the  discharge  first  half-period  (ti  =  8.9  ps)  exceeds  substantially  the  duration 
of  the  first  half-period  for  the  typical  APPT’s  with  the  oscillatory  discharge  (xi  ~  3  ps).  For 
the  APPT  under  study  maximum  current  value  Jmax  is  achieved  at  t=3.4  ps,  when  plasma  blob 
is  in  the  channel  zone  between  the  Teflon  bars  and  maximum  values  for  Ne  are  reached  in  the 
same  zone  also.  After  that,  when  plasma  comes  close  to  the  exit  plane  of  bars  (t  >  4  ps),  the 
discharge  current  drops  and  electron  density  decreases. 


215 


ISTC _ _ _ _ Partner  Protect  #2317v 

FINAL  TECHNICAL  REPORT 

J,  kA 


Figure  HI.  3  3  Time  dependence  of  discharge  current  for  APPT  model 
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Figure  111.34.  Time-dependent  distribution  of  maximum  electron  density  along 
channel  axis  of  the  next  generation  APPT 

With  subsequent  plasma  motion  through  the  interelectrode  gap  after  the  bar  exit  plane, 
electron  density  grows  during  the  essential  part  of  the  first  half-period  (t  >  4.8  p.s).  This  is 
evidently  connected  with  the  discharge  distribution  along  the  whole  electrode  length  and 
formation  of  outflowing  (“carry-out”)  currents  at  the  discharge  chaimel  outlet.  As  a  result, 
current  density  j  and,  correspondingly,  electron  density  Ne  grow  near  the  electrode  cut;  this 
agrees  with  the  data  of  j  measuring  by  magnetic  probes  presented  in  Electron  density 
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outside  the  channel  after  the  electrode  cut  drops  (Figure  111.34),  and  this  is  explained  by  the 
blob  expansion  after  it  has  left  APPT. 

Comparison  for  the  electron  density  time-dependent  distributions  in  APPT  of  next  generation 
and  typical  APPT  is  illustrated  by  Figure  111.35.  It  is  seen  that  in  the  APPT  of  next  generation 
and  typical  APPT  the  distributions  of  Ne  differ  from  each  other.  In  the  APPT  of  next 
generation  when  plasma  is  moving  from  the  cut  of  Teflon  bars  to  the  cut  of  electrodes  the 
electron  density  is  increasing.  This  is  coimected  with  discharge  expanding  at  the  end  of  the 
first  half-period  and  with  increasing  of  current  density  close  to  cut  of  electrodes. 

In  the  case  of  typical  APPT  with  oscillatory  discharge  and  low  level  of  performance,  the 
current  is  localized  near  the  bar  cut  during  nearly  the  entire  first  half-period  of  discharge  and 
is  not  carried  out  into  the  interelectrode  gap  within  the  discharge  charmel  part  from  the  bar  cut 
to  the  electrode  cut.  At  this  part  of  the  typical  APPT  discharge  channel,  electron  density  and 
current  density  drop  during  the  plasma  blob  motion  towards  the  APPT  outlet. 

Thus,  test  data  on  N©  measuring  witness  of  the  current  redistribution  in  the  channel  of  next 
generation  APPT.  With  this,  considerable  part  of  current  is  not  localized  in  the  area  of  bar  cut, 
but  is  carried  out  into  the  interelectrode  gap  after  the  bar  cut  providing  rather  efficient 
acceleration  of  propellant  by  electromagnetic  forces.  Such  discharge  dynamics  exactly  and 
such  distribution  of  electron  density  and  current  density  promotes  the  APPT  efficiency 
growth. 

Experiments  showed  that  current  sheet  propagates  inside  the  channel  at  an  average  velocity  of 
about  6-8  km/s.  Taking  into  account  that  the  velocity  registered  at  the  electrode  cut  is  over  20 
km/s  [1 1],  it  is  evidently  possible  to  consider  die  process  of  plasma  blob  acceleration  as  rather 
efficient. 

According  to  the  data  obtained  by  interferometric  method  in  the  area  after  the  Teflon  bar  cut, 
average  velocity  of  plasma  layer  with  maximum  value  of  Ne  near  the  anode  (~30  km/s)  is 
higher  substantially  than  that  near  the  cathode  (~20  km/s).  Thus,  inclination  of  lines 
coimecting  the  points  with  maximum  values  of  Ne  near  the  anode  and  cathode  is  observed  in 
the  discharge  channel,  and,  correspondingly,  the  inclination  of  electric  current  j  lines  relative 
to  the  electrodes.  This  distinction  in  the  distribution  of  Ne  and  j  is  evidently  coimected  with 
the  Hall  effect  exhibition. 

Figure  III.36  shows  the  Ne  variation  along  the  channel  length  during  the  plasma  movement 
near  the  cathode  (curve  1)  and  near  the  anode  (curve  2).  Distributions  of  Ne(x,  t)  are  not 
monotonous:  in  the  middle  part  of  the  zone  between  Teflon  bars  and  at  the  electrode  cut  the 
electron  density  is  higher  near  the  anode  than  near  the  cathode,  and  at  the  discharge  channel 
part  ranging  from  the  bar  cut  to  the  electrode  cut  the  values  of  electron  density  in  the  near¬ 
cathode  region  exceed  the  Ne  values  near  the  anode. 

The  indicated  differences  in  the  Ne  distributions  in  the  near-cathode  and  near-anode  areas  are 
evidently  connected  with  influence  of  competing  transversal  forces  -  electromagnetic  and 
gasdynamic  -  upon  plasma  in  the  discharge  channel.  Transversal  electromagnetic  force 
appears  as  a  result  of  Hall  effect  and  interaction  of  longitudinal  Hall  current  and  transversal 
magnetic  field  formed  by  the  APPT  operating  current. 
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Figure  IIL35.  Comparison  for  electron  density  time-dependent  distribution  along 
channel  axis  in  the  next  generation  APPT  and  the  typical  APPTl  -  the  next  generation 
APPT;  2  -  the  typical  APPT 
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Figure  IIL36.  Time-dependent  distributions  of  maximum  electron  density 
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Transversal  electromagnetic  force  presses  plasma  flow  to  the  cathode  and  is  counterbalanced 
by  the  transversal  gradient  of  pressure.  As  a  result,  gasdynamic  pressure  grows  near  the 
cathode,  and  plasma  blob  may  be  forced  out  from  the  cathode  and  deflected  to  anode  during 
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its  motion  under  the  influence  of  gasdynamic  force.  Thus,  dynamics  of  transversal  force 
influence  upon  plasma  in  the  discharge  channel  is  evidently  exhibited  by  the  periodic 
deflection  of  flow  to  cathode  and  anode.  It  should  be  noted  that  the  difference  in  Ne  values  in 
the  next  generation  APPT  channel  close  to  cathode  and  anode  near-cathode  and  near-anode 
areas  in  the  zone  between  bars  is  observed  in  the  case  of  time  intervals  of  about  1  ps  from  the 
discharge  start,  and  Hall  effect  starts  to  influence  substantially  the  dynamics  of  discharge 
evolution  at  the  very  early  stages  of  this  process  already. 

Plasma  blob  deflection  at  the  APPT  outlet  from  the  channel  axis  to  the  cathode  to  the  angle  of 
~  6°  under  the  influence  of  transversal  electromagnetic  force  is  coimected  with  the  Hall  effect 
exhibition  also.  Such  blob  deflection  was  discovered  while  measuring  the  flow  boundaries 
using  screens  [11]  and  recorded  while  fulfilling  this  work  also. 

Figure  III.37  shows  the  Ne  distribution  transverse  the  plasma  blob  at  a  distance  of  5  cm  from 
the  APPT  electrode  cut  (lines  IV-IV,  Figure  III.30).  Ne  distribution  is  non-symmetric  relative 
to  the  direction  of  plasma  blob  motion  (under  the  angle  of 6°  to  the  APPT  axis),  while  in  the 
blob  part  being  close  to  the  cathode  the  Ne  values  are  higher  than  near  the  anode. 
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Figure  111.37.  Distribution  of  maximum  electron  density  across  plasma  blob  at  the 
distance  of  5  cm  from  APPT  electrodes 


III.4.4  Piezoprobe  Measurements  of  Impact  Pressure  in  APPT  Plasma 

Flow 

The  main  purpose  of  the  piezoprobe  measurements  is  to  determine  the  "spread"  of  mass 
according  to  velocity  and  angle  influencing  the  operation  efficiency  of  APPT  as  a  thruster  and 
effecting  the  fields  with  variable  effectiveness  of  plasma  acceleration. 
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The  laboratory  model  of  the  next  generation  thruster  has  been  investigated.  Teflon  used  as  a 
plopellant.  The  energy  stored  in  the  capacitor  bank  was  W  =  100  J,  die  thrust  impulse  Pint  = 
2.3  mN'S,  the  pulse  mass  flow  rate  mpuUe  =  0.12  mg/pulse,  the  thrust  efficiency  iit  =  0.22  and 
the  specific  impulse  kp  =  19  km/s.  The  photo  of  thnister  is  presented  in  Figure  III.38.  The 
scheme  of  thruster  discharge  channel  is  presented  in  Figure  III.39. 

Piezosensor  of  impact  pressure  (piezoprobe)  is  a  thin  plate  of  piezoceramics  glued  between 
two  rods  made  of  quartz  glass.  The  length  of  the  rods  (100  mm)  has  been  chosen  on  the  basis 
of  the  fact  that  the  propagation  time  of  forward  and  reflected  acoustic  waves  in  quartz  was 
considerably  more  in  comparison  with  the  characteristic  time  of  the  plasma  flow  being  for  the 
given  thruster  type  t  ~  20  ^s.  Such  a  ratio  makes  possible  to  avoid  the  noise  effect  on  a 
recorded  signal  the  source  of  which  is  a  thruster  itself  and  its  systems.  Besides,  a  piezoprobe 
was  placed  into  a  copper  tube  to  provide  an  additional  noise  shielding;  and  a  coaxial  cable 
connecting  a  probe  wiA  a  pressure-seal  feed-through  of  a  vacuum  chamber  was  arranged  into 
a  flexible  screen,  a  kind  of  a  copper  braiding.  The  diameter  of  a  piezoprobe  end  part  is  4  mm 
to  sense  plasma  impact  pressure.  The  piezoprobe's  time  resolution  is  0.1  /is,  the  space 
resolution  is  5  mm  and  the  sensitivity  is  3-10*  Pa/V.  The  scheme  of  piezoprobe  is  presented  in 
Figure  III.40 

During  the  experiment  the  probe  was  moving  with  the  use  of  a  remote  controlled  coordinating 
device.  A  digital  two-beam  storage  oscilloscope  was  used  to  record  the  discharge  current 
curve  and  a  signal  from  a  piezoprobe.  One  of  the  oscillograms  is  presented  in 
Figure  iii  .41. 

The  signal  recorded  by  a  piezoprobe,  in  case  of  supersonic  plasma  flow,  is  proportional  to  the 
total  gas-dynamic  pressure  (p+pV^/2),  where  p  -  static  pressure;  p  and  V  -  flow  density  and 
velocity  respectively. 

The  total  signal  integral  in  any  plane,  perpendicular  to  a  flow  longtitudinal  axis  X, 
+  pV^/2)dydz  is  proportional  to  the  flow  impulse  in  this  plane  at  the  given  moment. 

The  integral  J  Pdt  according  to  the  whole  period  of  the  discharge  time  is  proportional  to  the 
total  thrust  impulse  Pint  created  by  the  thruster. 

The  measurements  were  performed  in  three  cross-sections,  corresponding  to  the  Teflon  bars 
cut  section  (1),  electrodes'  cut  section  (2)  and  distence  14  cm  from  the  electrodes'  cut  section 
(3).  The  positions  of  these  cross-sections  are  presented  in  Figure  III.42.  The  piezoprobe  was 
moving  along  the  axis  Y  and  Z  with  a  spacing  of  1  cm. 

The  result  of  measurements  in  three  cross-sections  are  presented  in  three-dimension  diagrams 
in  Figure  III.43-  Figure  111.64.  It  is  seen  that  in  all  cross-sections  the  maximum  flow  impulse 
is  shifted  to  the  thruster  anode. 

The  piezoprobe's  signal  integration  was  obtained  in  two  cross-sections'  plane  corresponding  to 
the  electrode's  cut  (X=0)  and  the  distance  14  cm  from  the  electrode's  cut  (X=14  cm).  The 
relevant  curves  P(0)  and  P(14  cm)  are  given  in  Figure  III.65.  There,  for  the  time  fixing  there 
is  a  discharge  current  J  oscillogram.  It  is  obvious  that  the  curve  P(0)  is  of  a  two-humped 
character.  The  second  maximum  of  the  curve  corresponding  to  a  flow  accelerated  within  the 
second  discharge  current  half  period  is  significantly  lower  by  3-4  times  than  the  first 
maximum  which  appropriates  to  the  first  half  period.  In  the  curve  P(14  cm)  there  is  no  second 
maximum.  It  can  be  explained  by  an  essential  larger  flow  divergence  according  to  the  angle  in 
the  second  half  period,  that  is  seen  in  high-speed  discharge  photograms.  The  curve  P(14  cm) 
is  also  much  smoother  than  the  curve  P(0),  that  is  explained  by  the  flow  inhomogeneity 
according  to  the  particles'  velocities. 

It  is  possible  to  determine  the  maximal  Vmax,  minimal  Vmin  and  middle  Vmid  particles' 
velocities: 

Vmax  =  14  cm  /  2  /IS  =  7-10^  m/s; 

Vmin  =  14  cm  /  8  /IS  =  1.75-10^  m/s; 
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Vmid  =  14  cm  /  4  j[is  =  3,5-10'*  m/s. 

The  average  angle  of  flow  divergence  (only  in  die  first  half  period)  can  be  estimated  by  the 
relation  of  integral  impulses  J  Pdt  in  cross-sections  (X=0)  and  (X=14  cm).  The  appropriate 

evaluations  give  angle  2(p  ~  28°. 


Figure  III.38.  APPT  laboratory  thruster 


221 


ISTC 


Partner  Pro  ject  it2317p 


FINAL  TECHNICAL  REPORT 
Z 


Figure  III.39.  APPT  discharge  channel 
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Figure  ni.40.  Piezoprobe 
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Figure  III.41.  APPT  discharge  current  J  and  piezoprobe's  signal  P  curves 
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Figure  III.42.  Piezoprobe  measurement's  cross-sections 
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Figure  111.43.  Distribution  of  impact  pressure  in  Teflon  bars  cut  section  (cross-section  1) 
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Figure  ni.44.  Distribution  of  impuct  pressure  in  Teflon  bars  cut  section  (cross-section  1) 
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Figure  III.45.  Distribution  of  impact  pressure  in  Teflon  cut  section  (cross-section  1) 


Figure  in.46.  distribution  of  impact  pressure  in  Teflon  bars  cut  section  (cross-section  1) 
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Figure  III.48.  Distribution  of  impact  pressure  in  Teflon  cut  section  (cross-section  1) 
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Figure  III.49.  Distribution  of  impact  pressure  in  Teflon  cut  section  (cross-section  1) 


Figure  III.50.  Distribution  of  impact  pressure  in  Teflon  cut  section  (cross-section  1) 
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Figure  III.51.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 


Figure  111.52.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  ni.53.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  111.54.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  111.55.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  III.56.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  III.57.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  III.58.  Distribution  of  impact  pressure  in  electrodes'  cut  section  (cross-section  2) 
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Figure  ni.59.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  III.60.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  III.61.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  III.62.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  ni.63.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  111.64.  Distribution  of  impact  pressure  at  a  distance  14  cm  from  the  electrodes'  cut 

section  (cross-section  3) 
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Figure  III.65.  Integral  piezoprobe's  signal  P  in  cross-sections  2  (0  cm)  and  3  (14  cm) 


III.4.5  Ultra  high  speed  photography  with  visualization  of  plasma  flow 

Ultra  high-speed  photography  of  discharge  and  experimental  determination  of  Mach  numbers 
in  plasma  flow  were  tried  out  using  the  measurement  procedure  for  the  APPT  of  next 
generation.  The  appearance  of  thruster  laboratory  model  is  presented  in  Figure  III.38.  For 
increasing  the  discharge  glowing  the  energy  stored  in  the  capacitor  bank  was  increased  from 
W  =  100  J  to  W  =  150  J.  In  this  case  die  APPT  thrust  impulse  was  Pint  =  4.5  mN«s,  the  pulse 
mass  flow  rate  mpuise  =  0.22  mg/pulse,  the  flirust  efficiency  rjt  =  0.30  and  the  specific  impulse 
was  Isp  =  19  km/s.  The  oscillogram  of  discharge  current  is  presented  in  Figure  111.66.  High¬ 
speed  photographic  record  of  plasma  flow  was  performed  with  the  help  of  ultra  high-speed 
photographic  camera.  The  camera  has  a  long-focus  objective,  rotating  mirror  and  an 
electromagnetic  or  electrodynamic  shutter.  The  maximum  speed  of  mirror  rotation  is  60,000 
rev/min.  The  temporal  resolution  is  ~2  fjis;  the  spatial  resolution  is  ~1.5  mm;  the  time  spacing 
between  frames  is  ~2  /is. 

To  make  the  direction  and  velocity  of  plasma  flow  visible  a  sharpened  wedge  of  quartz  glass 
was  put  into  a  flow.  The  wedge  thikness  was  1  mm,  the  width  -  20  mm,  the  angle  of  point 
sharpening  ~30®.  The  wedge  was  moving  in  the  flow  according  to  two  axis  (X  and  Y)  with  the 
help  of  a  remote  controlled  coordinate  device.  In  this  case  the  shocks  typical  for  a  supersonic 
flow  were  generated  at  the  wedge  sharp  point  and  were  easily  observed  in  photograms.  The 
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vector  direction  of  the  flow  speed  and  Mach  number  M=l/sin(a/2),  where  a  is  an  angle 
between  oblique  shocks,  can  he  determined  by  the  shock  slope  to  a  longitudinal  axis  X. 

The  photographic  flow  record  was  carried  out  according  to  two  mutually  perpendicular 
optical  axis  Y  and  Z,  perpendicular  to  longitudinal  thruster  axis  X.  The  wedge  was  positioned 
by  die  coordinate  device  in  four  cross-sections,  corresponding  to  the  Teflon  bars  cut  section 
(1),  middle  of  electrodes'  length  (2),  electrodes'  cut  section  (3)  and  distance  7  cm  from  the 
electrodes'  cut  section  (4).  The  positions  of  these  cross-sections  are  presented  in  Figure  III.67. 
A  characteristic  pictures  of  oblique  shocks  is  presented  in  photographs  in  Figure  III.68,  which 
was  taken  by  a  mirror  photographic  camera  in  the  integral  discharge  glowing. 

Ultra  high-speed  photograms  of  thruster  discharge  and  shocks  in  plasma  flow  are  presented  in 
Figure  111.69  -  Figure  III.78.  The  bright  glowing  spot  visible  in  photograms  below  the 
discharge  channel  is  a  beam  of  the  laser  used  for  high-speed  camera  optical  tuning. 

In  photograms  in  Figure  111.69  -  Figure  III.76  the  optical  axis  of  the  objective  was  directed 
along  the  axis  Y.  It  is  seen  the  discharge  expansion  from  the  Teflon  bars  cut  section  (cross- 
section  1)  to  the  electrodes'  cut  section  (cross-section  3).  The  discharge  glowing  is  beginning 
near  anode  in  t  ~  1  /is  and  finishing  in  t  ~  15  ^s.  The  maximum  of  glowing  is  registering  in  t  ~ 
5  jUS,  that  is  about  2  /is  later  than  the  maximum  of  discharge  current.  The  oblique  shocks  near 
the  wedge  sharpening  are  visible  from  t  ~  3  /is  to  t  ~  13  /is. 

In  photograms  in  Figure  in.77  -  Figure  ni.78  the  optical  axis  of  the  objective  was  directed 
along  the  axis  Z.  In  photograms  in  Figure  III.77  it  is  seen  that  a  flow  has  two  characteristic 
glowing  areas: 

1)  the  main  flow  core,  directed  along  the  longtitudinal  axis  with  low  divergence, 

2)  a  peripheral  flow  directed  along  the  end  surfaces  of  the  Teflon  bars  with  a  large  divergence 
angle. 

The  flow  core  is  inclined  to  the  longitudinal  axis  at  the  angle  of  P  ~  5®  in  the  cathode 
direction.  In  the  first  discharge  half  period  the  flow  core  is  flie  main  source  of  glowing.  In  the 
second  half-period  the  core  is  split  up  and  the  peripheral  flow  becomes  the  main  source  of 
glowing. 

The  pictures  of  oblique  shocks  obtained  in  photogram  in  different  flow  points  and  different 
moments  of  time  are  given  in  Figure  III.79  -  Figure  III.90.  It  is  seen  that  the  Mach  number 
near  anode  is  essential  less  (M  =  1.5  -  1.8)  than  the  Mach  number  in  another  parts  of  the 
plasma  flow.  In  the  flow  core  the  Mach  number  M  =  2.5  ±  0.2  practically  constant  within  the 
first  discharge  half  period  is  a  distinctive  feature  of  the  plasma  flow  created  by  APPT. 

The  calculated  estimations  testify  that  this  number  M  corresponds  to  the  ratio  of  flow  velocity 
V  to  the  AlfVen  velocity:  Ma  =  V/Ca,  rather  than  the  ratio  of  flow  velocity  V  to  the  thermal 
sonic  velocity:  Mt  =  V/ct  (it  was  calculated  that  in  APPT  plasma  flow  Mt  >  10).  It  means  that 
in  flie  first  half  period  plasma  has  a  rather  high  degree  of  ionization,  as  a  result  of  it,  Alfven 
magnetosonic  waves  are  considered  to  be  the  basic  mechanism  to  transfer  the  disturbances 
instead  of  thermal  sound. 

We  can  estimate  flie  angle  of  flow  divergence  in  planes  XY  and  XZ  by  a  bisectrix  inclination 
angle  of  oblique  shocks  to  the  thruster  longitudinal  axis  X.  The  corresponding  measurements 
in  photograms  showed  that  in  the  first  discharge  half  period  this  angle  did  not  exceed  2(p  ~ 
20®.  In  the  second  half  period  plasma  glowing  is  not  enough  to  make  photographic  record  of 
shocks. 
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Figure  III.68.  APPT  discharge  and  oblique  shocks  integral 
photo  (a  -  cross-section  3,  b  -  cross-section  4) 
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Figure  III.69.  Photogram  of  APPT  discharge  and  shocks  (optical  axis 
Y,  cross-section  1) 
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V,  cm  t  =  3  mcs 


Figure  III.79.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 


y,  cm  t  =  5  mcs 


Figure  111.80.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 
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Figure  III.81.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 
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Figure  III.82.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 
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y,  cm  t  =  11  mcs 


Figure  III.83.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 


y,  cm  t  =  13  mcs 


Figure  III.  84.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Y) 
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Figure  111.85.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Z) 


y,  cm  t  =  5  mcs 


Figure  III.86.  Diagram  of  oblique  shocks  in  APPT’s  plasma  flow  (axis  Z) 
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Figure  III.87.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Z) 


y,  cm  t  =  9  mos 
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Figure  111.88.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Z) 
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Figure  III.89.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Z) 
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Figure  IIL90.  Diagram  of  oblique  shocks  in  APPT's  plasma  flow  (axis  Z) 
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III.5  Development  and  tests  of  the  APPT  model 


For  a  number  of  years  after  the  effect  of  rail  APPT  thrust  efficiency  substantial  growth  under 
certain  conditions  had  been  revealed,  main  attention  was  focused  on  the  study  and 
experimental  refinement  for  the  models  of  highly  efficient  thruster  operating  within  the  power 
consumption  range  of  (60-150)  J.  Such  approach  seemed  to  be  justified,  because  to  our  mind 
within  ftis  power  range  exactly  the  APPT  is  the  most  universal  thruster  capable  to  solve 
problems  of  attitude  control  for  a  small  satellite  of  50-500  kg  in  mass  most  effectively.  A 
number  of  reports  devoted  to  this  subject  were  made  at  international  conferences  and  several 
articles  concerning  this  issue  were  published.  Some  results  of  investigations  and  refinements 
for  the  models  with  the  stored  energy  of  100  J  and  150  J  with  higher  efficiency  than  that  of 
traditional  thrusters  are  presented  in  the  above  parts  of  the  present  report  also. 

An  opinion  was  already  expressed  at  the  initial  stage  of  highly  efficient  APPT  study,  that 
without  the  introduction  of  additional  matching  devices  (it  is  extremely  undesirable  from  the 
point  of  view  of  securing  the  thruster  competitiveness),  the  external  electric  circuit  of  real 
APPT  have  a  limited  abilities  for  its  optimization.  In  particular,  estimations  showed  that  when 
the  energy  stored  in  the  accumulator  is  decreased  down  to  (20-30)  J,  potentialities  for 
improving  the  thruster  performance  by  matching  the  elements  of  its  circuit  might  appear  to  be 
exhausted. 

In  view  of  the  current  vast  interest  to  APPT  with  the  stored  energy  of  up  to  20  J,  that  is 
capable  to  provide  efficient  orbital  and  attitude  control  for  micro-  and  nano-satellites,  we 
made  an  attempt  for  experimental  determination  of  the  minimum  stored  energy,  at  which  the 
thruster  efficiency  could  be  increased  by  optimizing  the  elements  of  its  circuit.  For  this,  APPT 
laboratory  models  with  the  nominal  energy  of  40  J  and  30  J  were  developed,  which  allowed 
the  energy  variation  within  25%  of  the  nominal  value  by  changing  the  discharge  voltage. 
Besides,  at  the  discoimection  of  one  of  three  capacitors,  the  30  J  model  could  be  transformed 
easily  into  the  20  J  model.  Thus,  the  range  of  stored  energy  from  50  J  down  to  15  J  could  be 
covered  during  ftie  test.  External  appearance  of  ftie  models  is  presented  in  Figure  111.91  and 
Figure  111.92,  correspondingly.  The  model  with  the  stored  energy  of  40  J  (Figure  III.91)  is 
placed  on  a  thrust-metering  device  designed  specially  for  measuring  low  thrust.  Energy 
storage  of  both  models  comprises  several  capacitors,  while  the  discharge  chaimel,  formed  by 
two  plane  electrodes  and  propellant  bars  fed  from  the  side  directions,  is  located  at  a  minimum 
distance  from  the  storage. 

Results  of  experimental  refinement  for  APPT  models  are  presented  in  the  Table  4. 
Performance  data  for  the  APPT  models  with  the  standard  power  supply  circuit  is  from 
mainly.  It  is  obvious  from  data  presented  in  Table  4,  that  within  the  energy  range  from  50  J 
down  to  20  J  at  least  the  APPT  with  improved  power  supply  circuit  has  considerably  better 
characteristics  than  the  thruster  with  the  standard  circuit.  At  that,  within  the  above  energy 
range  the  relation  of  values,  for  the  thruster  efficiency  for  example,  which  is  close  to  two,  is 
kept  practically  always.  At  the  energy  value  of  below  20  J  only,  some  convergence  of 
efficiency  values  may  take  place. 
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Figure  III.91.  APPT  40  J  appearance 


Figure  III.92.  APPT  20  J  appearance 


We  did  not  manage  to  find  traditional  analogue  (prototype)  for  the  improved  15  J  thruster  and 
to  compare  their  characteristics,  but  test  results  obtained  give  good  reason  to  believe  that 
advances  of  highly  efficient  thruster  will  be  kept  down  to  energy  level  of  15  J  at  least. 
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Table  4 


Storage 

energy 

Thrusters  with  improved  circuit 

Thrusters  with  standard  circuit 

Impulse 

bit 

Propellant 

consumption 

Thrust 

efficiency 

Propellant 

consumption 

Thrust 

efficiency 

J 

mNs 

kg 

% 

kg 

% 

50 

1.25 

0.8  -lO’’ 

0.2 

0.85 

0.65  -lO'’ 

0.11 

40 

1.0 

0.7  •10-'^ 

0.18 

0.74 

0.62-10'’ 

0.1 

30 

0.7 

0.6  lO'^ 

0.13 

0.46 

0.4-10'’ 

0.08 

20 

0.47 

0.6  -lO’’ 

0.11 

0.3 

0.35- lO’’^ 

0.06 

15 

0.35 

0.4  lO'’ 

0.09 

no  data 

It  follows  from  above  presented  that  investigations  of  highly  efficient  APPT  within  the  range 
of  relatively  low  energies  are  of  preliminary  nature  and  should  be  continued  for  obtaining  the 
array  of  test  points  sufficient  for  making  more  definite  conclusions.  In  particular,  in  the 
nearest  future  it  is  planned  to  continue  experimental  comparison  for  the  traditional  and 
improved  thruster  models  within  the  range  of  relatively  low  energies.  APPT  models  with  the 
stored  energy  of  15  J  and  10  J  of  both  designs  will  be  produced  and  tested  for  this. 
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III.6  Conclusion 

Two-dimensional  semi-phenomenological  numerical  model  of  APPT  have  been  developed. 
Boundary  condition  for  propellant  flow  rate  is  accepted  from  experiments  and  calculations  on 
tire  base  of  propellant  thermal  degradation.  Developed  model  adequately  describes  formation 
and  acceleration  of  a  plasma  flow  in  coaxial  accelerating  channel.  The  simulated  integrated 
parameters  of  considered  APPT  have  appeared  close  to,  measured  in  experiments.  The  model 
gives  a  correct  qualitative  picture  of  separate  stages  of  discharge  in  APPT:  plasma  flow 
formation,  appearance  of  tubular  stmcture  at  the  end  of  the  central  electrode,  focusing  of  a 
flow  etc.  The  features  of  propellant  flow  rate  determine  monotonous  character  of  a 
distribution  of  an  electron  density  and  especially  an  axial  velocity  of  a  plasma,  and  magnetic 
field  in  accelerating  channel.  The  developed  numerical  model  of  APPT  can  serve  as  a  tool  to 
improve  thruster  parameters  through  experimentation  -calculation  process. 

Energy  flux  onto  propellant.  Teflon  mass  loss,  current  and  voltage,  near-propellant  bar 
electron  density  have  been  measured  in  APPT  stand.  Energy  fraction  released  from  the 
discharge  region  upon  the  propellant  surface  in  the  APPT  is  near  10'^  of  the  energy  stored  in 
the  power  supply  source.  Main  part  of  measured  flux  is  transferred  by  particles.  For  bank 
energy  20-  100  J,  measured  energy  in  a  thermal  skin  layer  is  less  than  accounted  from  model 
of  Teflon  thermal  degradation.  There  is  no  observed  delay  between  the  current  beginning  and 
plasma  appearance  near  propellant.  So,  at  the  beginning  of  the  discharge,  propellant  flow  rate 
is  provided  due  to  non-thermal  or  high  power  energy  flux  mechanisms.  Near  propellant  the 
electron  density  depends  on  the  input  voltage  and  measured  at  the  level  of  lO’^cm'^  for  0.5  kV 
(inter-electrode  gap  voltage)  with  some  deviation  due  to  shot-to-shot  changes  in  thruster 
firing.  Thickness  of  boundary  layer  is  order  of  0.1  mm 

Electron  density  measurements  in  the  APPT  discharge  chaimel  in  the  area  between  propellant 
(Teflon)  bars  and  in  the  interelectrode  gap,  as  well  as  outside  the  discharge  channel,  were 
made  by  the  method  of  laser  interferometry.  A  diagnostic  system  developed  on  the  basis  of 
three-mirror  laser  interferometer  was  used  during  tests.  Due  to  the  optimization  of  the  electric 
circuit  parameters,  the  studied  APPT  has  high  specific  characteristics.  In  particular,  current 
efficiency  and  specific  impulse  at  the  stored  energy  of  100  J  comprise  24%  and  18km/s, 
correspondingly. 

Studies  for  the  spatial  and  temporal  distributions  of  the  electron  density  in  the  discharge 
chaimel  make  possible  to  reveal  main  features  of  physical  processes  of  the  next  generation 
APPT’s  possessing  the  improved  thruster  characteristics.  Dynamics  and  structure  of  plasma 
flow  in  high-efficiency  APPT  with  stored  energy  of  150  J  were  studied. 

The  data  confirming  the  established  concept  in  respect  to  the  dominant  role  of 
electromagnetic  mechanism  under  plasma  acceleration  in  high-efficiency  APPT  were 
obtained  by  using  the  ultra-high  speed  photo  camera.  These  measurements  have  shown  that 
plasma  flow  in  such  thmster  have  the  insignificant  divergence  near  the  electrode  cut.  The 
insignificant  part  of  propellant  flies  apart  at  the  more  angle  relevant  the  thermal  velocity 
under  the  discharge  decay  only.  This  propellant  part  can  constitute  some  threat  to  optical 
surfaces  of  spacecrafts. 

Piezoprobe  measurements  of  plasma  flow  impulse  distributions  have  shown  that  the  impulse 
maximum  is  displayed  to  the  thruster  anode.  This  must  be  taken  into  account  in  the  design  of 
propulsion  systems  based  on  APPT.  The  divergence  angle  of  plasma  flow,  which  was 
obtained  at  piezoprobe  measurement  processing,  was  somewhat  higher  (about  30°)  in 
comparison  with  measurement  carried  out  by  ultra-high  speed  photo  camera.  However,  in  this 
case  the  good  focusing  of  flow  was  observed. 
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Two  APPT  models  having  the  level  of  bank  energies  20  J  and  40  J,  are  developed  and  tested. 
The  improved  APPT  characteristics  are  obtained  with  completion  of  an  optimization  between 
parameters  of  an  electric  circuit  and  the  sizes  of  an  accelerating  channel  of  the  thruster,  when 
the  discharge  is  close  to  aperiodic.  Laboratory  refinement  of  these  models  have  shown  that 
engineering  approach  making  possible  to  increase  the  APPT  operation  efficiency  may  be  used 
over  the  energy  range,  which  is  limited  on  the  underside  by  energy  values  close  to  20  J. 

Based  on  interest  which  is  aroused  by  APPT  with  stored  energy  of  ~20  J  and  less  it  would  be 
useful  to  conduct  extended  studies  for  the  purpose  of  raising  the  thruster  operation  efficiency 
over  this  energy  region. 

Two  presentations  have  been  done  under  Project  in  International  Electric  Propulsion 
conference,  IEPC-2003,  Toulouse,  France: 
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III.7  Annex  1.  Propellant  Flow  Rate  Modeling 

III.7.1  Introduction 

Almost  all  APPT's  to  date  have  used  solid  Teflon  (Polytetrafluoroethylene,  or  PIPE)  as  a 
solid  propellant,  because  it  does  not  exhibit  carbon  buildup  on  itself  or  the  surrounding 
structure  and  because  its  high  temperature  degradation.  In  connection  with  this,  the  dynamics 
of  the  insulator  ablation  is  of  a  definite  interest  for  comprehending  the  plasma  acceleration 
mechanism  in  similar  instruments  and  for  finding  an  optimal  mode  of  their  operation. 

In  the  APPT,  the  considerable  energy  fraction  released  from  the  propellant  surface  may  be 
transferred  not  only  by  particles  but  ultraviolet  radiation  in  the  wavelength  range  near  1000  A. 
The  depth  of  this  radiation  absorption  in  polymers  does  not  exceed  10'^  cm;  meanwhile  the 
thickness  of  an  propellant  material  heating  due  to  the  heat  conduction  is  about  10"^  cm.  At 
present  the  relative  role  of  conductive  and  irradiative  mechanisms  for  energy  transfer  to  the 
propellant  is  not  exactly  determined.  Now  conductive  mechanism  is  more  expanded  and 
developed.  Probably  the  main  mechanism  of  energy  transfer  to  propellant  depends  on  thruster 
operation  mode.  In  any  case,  if  one  assumes  that  the  energy  flux  released  upon  the  propellant 
surface  from  the  discharge  zone  and  that  the  evaporated  substance  is  instantaneously 
removed,  the  problem  of  propellant  ablation  in  the  APPT  will  be  reduced  to  the  solution  of  a 
nonlinear  heat  conduction  equation  with  a  moving  interface.  In  this  connection,  studying  the 
propellant  ablation  in  these  instruments,  one  is  forced  to  consider  die  kinetic  of  thermal 
degradation  that  considerably  complicates  the  problem.  On  the  other  hand,  the  thickness  of  an 
evaporating  layer  in  the  propellant  is  very  small,  and  therefore  one  can  limit  himself  by  a 
study  of  ID-heat  conduction  equation. 

Below  we  take  into  account  thermal  PTFE-degradation  kinetics  to  solve  problem  of  propellant 
ablation  in  the  APPT.  The  approach  of  the  PTFE-  degradation  kinetics  developed  by 
Madorski  is  expounded.  The  statement  and  the  results  of  solving  the  problem  are  given, 
taking  account  of  polymer  degradation  obeying  the  kinetic  equation  of  the  first  order.  Also  it 
is  proposed  to  use  two  branches  of  degradation  reaction  more  to  represent  the  thermal 
polymer  degradation  at  the  energy  flux  densities  to  the  propellant  >10'*  W/cm^. 


I1I.7.2  Kinetics  of  Teflon  Thermal  Degradation 

Let  us  briefly  speak  about  modem  ideas  concerning  the  thermal  degradation  of  polymers. 
Since  the  propellants  made  of  Teflon  have  been  more  often  studied  in  the  APPT,  all  the 
quantitative  data  are  referred  to  this  polymer. 

According  to  the  following  mechanism  of  degradation  is  applicable  to  all  the  polymers 
which  undergo  the  thermal  degradation  up  to  a  monomer,  dimer  and  to  die  other  oligomeric 
molecules  produced  as  a  result  of  rupture  in  the  main  chain: 


Reaction 

Reaction  rate 

Energy 

Initiation  (formation  of  radicals)  Qn  ->  Ri  +  Rn-i 

2  ki  n  Qn 

ei 

Growth  of  a  chain  (product  formation) 

Rn  Rn-l+M 

k2R 

S2 

Rupture  of  the  chain  (recombination) 

RiH-  Rj— >Qi+i 

2k4R^ 

U  1 

Here,  Qn  is  the  concentration  of  polymer  molecules  including  n-  links  of  C2F4;  Ri  is  the 
concentration  of  radicals  having  i-  monomeric  links;  M  is  the  concentration  of  monomers. 
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If,  along  with  heating,  the  polymer  is  irradiated,  there  will  be  one  opportunity  more  to  initiate 
the  reaction 

hv 

Qn  — >R.+Rn-i 

The  rate  of  this  reaction  is  equal  2(pInQn  ,  where  I  is  die  radiating  power  released  under 
absorption  in  a  polymer,  (p  is  the  radiation-chemical  yield  of  radicals. 

As  for  the  radiation-chemical  yield  of  irradiation,  cp,  the  mechanism  of  electron  excitation  and 
of  the  subsequent  dissociation  results  in  that  the  radiation-chemical  yield  for  the  y  -radiation 
and  for  ultraviolet  one  turn  out  to  be  the  same  by  the  order  of  magnitude.  Taking  account  of 
the  fact  that  a  considerable  fraction  of  the  energy  released  upon  the  propellant  is  transferred 
by  radiation,  one  can  expect  that  the  contribution  of  radiation  into  the  reaction  of  initiation 
will  be  essential,  at  the  initial  stages  of  a  discharge  in  particular,  when  the  propellant 
temperature  is  not  great  and  radiation  initiation  of  the  depolymerization  reaction  can 
considerably  exceed  the  thermal  one. 

Thus,  taking  account  of  the  radiation  initiation,  the  equation  of  radical  formation  can  be 
written  in  the  form 

^=2(k,  +  ^i)pQ„-2k,R"  (2.1) 

On  the  other  hand,  according  the  chain  growth  reaction,  the  rate  of  monomer  formation  is 
equal 


dM  .  „ 
dt 


(2.2) 


Let  us  assume  that  the  molar  fractions  of  radicals  and  monomers  in  the  sample  are  small,  i.e. 

XiRi«Z!r>Qn  and  M  «XnQn.  (C«0). 


Making  transition  to  dimensionless  variables  -  relative  concentration  of  radicals  R->^R/(no(l- 
C)  and  conversion,  C=M/no,  one  obtains  a  set  of  equations  representing  the  thermal  Teflon- 
degradation  from  (2.1)  and  from  (2.2): 


dt 


=  2(k,  +  ^l)-2k4R"no(l-C), 


1  dC 
1-C  dt 


=k2R 


(2.3) 


And,  finally,  if  one  assumes  that  the  radiation  initiation  can  be  neglected  in  comparison  with 
the  thermal  one,  the  following  equation  of  the  pol)mier  degradation  will  take  place  for  the 
equilibrium  case: 
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(i-cy^  dt 

R  =(r^)' 

k4no 


'Mo 


0.5 


(2.4) 


The  equation  (2.4)  coincides  with  the  conversion  equation  produced  in  and  the  quantity  k2 
(ki  /  k4no)  is  a  constant  of  the  reaction  k"  =4.14  lO'*  exp(-3.5/Tev). 

Some  thermo-physical  parameters  and  kinetic  constants  for  the  Teflon  are  given  Table  1. 


Table.  1 .  Teflon  thermo-physical  parameters 


Density  at  a  normal  temperature,  p 

2.2  -  2.3  g/cm" 

Heat  capacity  at  normal  temperature,  Cp 

1.06  J/gradC 

Heat  conductivity  at  normal  temperature,  ti 

2.5  10'^  J/(cm  s) 

Energy  of  bonds: 

C-C 

2.78  eV 

c=c 

4.8  eV 

C-F 

5.0  eV 

Constants: 

thermal  initiation,  kj 

-2.110^*  exp(-^)  s'* 

'  ev 

chain  growth,  k2 

1 

710*’ exp(-  )s' 

'  ev 

chain  rupture,  k4 

1  3 

510'®  exp(--f-)  s'* 

'  ev 

radiation-chemical  yield  (order  of  magnitude)  (p 

io'-^  ey-'* 

III.7.3  Equation  of  Energy  Transfer  in  the  Propellant 

The  equation  of  energy  transfer  in  the  propellant,  in  the  immobile  coordinate  systems,  has  the 
form: 


ds  d  <?T 


(3.1) 


Here  e  is  the  intrinsic  energy  in  1  cm^  of  the  substance.  Taking  account  of  the  polymer 
degradation  mechanism  given  in  Section  2,  the  intrinsic  energy  of  the  substance  can  be 
written  in  the  form 


£  =  yo(CpT  +^Rno+ff2Cno), 

(3.2) 

* 

s  =  P  (CpT  -i--^R+£2C), 

(3.2a) 

where  =  4900  J/g  and  £2  =  1517  J/g. 
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Substituting  (3.2a)  into  (3.1),  one  obtains 


R+£2C)  = 


(3.3) 


Adding  the  relationships  (2.3),  (2.4)  in  the  equation  (3.3)  and  taking  account  of  the  fact  that 
C  »  0,  one  obtains  the  set  of  the  equations  which  can  be  used  for  representation  of  the  pulsed 
propellant  depol)mierization  (evaporation)  in  the  APPT 


d  s\  d  dJ 

P  2  ^  '  dY'  dY 

■^  =  2k,-2k4R2no,-^  =  k2R. 


(3.4) 


Before  formulating  the  boundary  conditions  for  (3.4),  let  us  make  some  assumptions  on  a 
polymer  degradation  model  in  oiu:  case.  First,  let  us  assume  that  the  density,  heat  conductivity 
and  heat  capacity  of  the  pol)mier  remain  to  be  constant.  Thus  the  polymer  degradation 
problem  is  reduced  to  the  problem  similar  to  the  evaporation  of  a  fluid  from  the  surface.  Since 
the  heated  layer  thickness,  5,  is  small  and  the  thickness  of  the  layer,  where  an  intense  polymer 
degradation  occurs,  is  ~  8/10,  one  can  assume  that  the  diffusion  of  a  monomer  from  the 
propellant  goes  on  at  the  infinite  rate.  At  high  depths  of  heating  the  monomer  yield  can  be 
controlled  by  diffusion  process  Second,  let  us  assume  that  the  maximal  propellant  surface 
temperature  does  not  exceed  1300°K,  i.e.  a  relative  amount  of  radicals  in  the  layer  is  small. 

The  relative  concentration  of  radicals  in  the  heated  layer  will  attain  its  equilibrium 
value  for  the  time,  t  ~  1/ki ,  which  is  about  2.  10'*  s  for  T=1200‘’  K  fiiat  is  considerably  lower 
than  the  characteristic  plasma  acceleration  time  in  the  instruments  under  consideration. 
Therefore  the  number  of  radicals  in  1  cm*  of  the  substance  is  determined  by  the  second 
relationship  in  (2.4),  and  the  contribution  introduced  by  radicals  into  the  intrinsic  energy,  e  ,  is 

considerably  smaller  than  Cp  T.  Indeed,  for  T  =  1200°  K,  CpT  =960  J/g  and  s*  —  =  105  - 

•  ^ 

125  J/g 

Taking  the  above  expressed  assumption,  the  set  (3.4)  can  be  written  in  the  form 


£I 

dO 

d\ 


—  k2  ( 


dY^ 
ki 


k4no 


e\  dQ> 

.)0.5  =4  4410’*  exp 

I  CM 


(3.5) 


Eqs.  (3.5)  are  solved  in  the  rage  s(t)  J  x  J  f ,  where  s  (t)  is  the  phase  interface.  The  boundary 
and  initial  conditions  for  (3.5)  have  the  form 


atx  =  s(t) 

T  (t,  00)  =  To , 
T(0,x)  =  To. 


(3.6) 
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The  interface  rate, 


can  be  found  from  the  relationship 


d\ 


ki 


k2(r-^r  dx. 


k4no 


(3.7) 


Eqs.  (3.5,  together  with  (3.6)  and  (3.7),  were  solved  numerically.  The  equations  were  made 
dimensionless..  The  following  parameters  DT  =  900®,  to  =  hCpT  DT^/Ho^  were  chosen  to  be 
scaling  law  quantities.  Having  remained  the  previous  designations  of  the  quantities,  one 
obtains  a  dimensionless  heat  conduction  equation 
l-A^C 


d\ 

where  A= 


Cp  AT 


Cp  AT  +  s\ 


A  di  ’ 

.  Initial  and  boundary  conditions,  in  this  case,  acquire  the  form 


dl 

dy. 


=  -l,  T(0,x)  =  a,  T(t,oo)  =  a, 


(3.9) 


where  a  =  AT-l/To.  The  monomer  3deld  rate  was  determined  from  the  relationship 
^  =  Bexp(-Y),  (3.10) 


where  B=  4.14  10’®to,  P  =45.1.  The  heat  conduction  equation  was  solved  numerically. 


III.7.4  High  Temperature  Branches  of  Degradation 

As  noted  above,  at  the  energy  flux  densities  upon  the  propellant  surface  greater  than 
5.10^J/cm^  s,  the  stationary  propellant  surface  temperatures  exceed  1250  K.  A  relative 
amount  of  radicals  in  the  substance  becomes  equal  to  ~  1  that  contradicts  to  the  assumption 
about  the  smallness  of  R,  made  in  Section  2.  Unfortunately,  there  are  no  experimental  data  on 
the  Teflon  -degradation  rate  at  the  temperatures  >  800  K. 

At  the  transition  to  higher  energy  flux  densities  upon  the  propellant  surface  and,  as  a  result,  to 
higher  polymer  degradation  temperatures  one  should  take  into  account  of  the  fact  that  the 
formation  of  difluoro-methylene  radicals  (CF2)  and  excited  molecules  (C2F4)  ,  at  such 
temperatures,  can  occur  with  greater  probability  than  the  formation  of  monomers  at  the  ends 
of  a  polymer  chain.  In  this  case,  the  chain  reaction  diagram  should  include  two  reactions 
more: 


Reaction 

Reaction  rate 

Formation  of  a  difluoro-methylene  radical 

Ri  ->  Ri-1/2  +  Ri/2 

kiR/2 

Formation  of  the  excited  molecule  of  monomer 

Ri-^Ri.i+M* 

kiR 

If  one  assumes  that  the  thermal  disrupture  of  bonds  of  the  ends  of  a  chain  occurs  with  the 
same  probability  as  that  at  its  other  pieces,  one  will  be  able  to  use  the  initiation  constant,  ki , 
for  the  estimations. 

The  Teflon  -  degradation  rate,  taking  account  of  all  the  above  -  considered  possible  reactions 
considerably  exceeds  the  monomer  formation  rate  on  radicals  at  the  temperatures  >  1 100  K. 
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At  the  same  time,  the  diagram  under  consideration  does  not  contradict  the  known 
experimental  results  since  die  contribution  of  these  two  chain  reaction  branches  is 
negligible  in  the  temperature  range  up  to  800  K. 

Let  us  assume  that,  in  this  case  under  consideration,  a  relative  amount  of  radicals  in  a 
substance  is  small.  Than,  the  heat  conduction  equation  will  acquire  the  form 

^  =  -TT  -l/Cr[3/2  ki8*i  +  k2S*2  ]  (ki/  k4no)°'^  (4.1) 

at  ox 

(values  of  the  constants  kl  ,k2  ,  k4  are  given  in  Section  2).  The  boundary  and  initial 
conditions  of  problem,  as  well  as  an  expression  for  the  interface  velocity,  remain  to  be 


previous. 

It  is  easy  to  obtain  that  the  temperature,  Tst ,  weakly  depends  on  Ho.  When  the  energy  flux  is 
changed  from  10^  to  10®  J/cm^  s,  the  temperature  varied  by  150°  only.  If  one  takes  account  of 
the  fact  the  degradation  mechanisms  introduced  in  this  Section  have  the  reaction  energy  Si , 
the  average  Teflon  -degradation  energy,  e,  will  be  determined  by  the  following  expression 

1.5 k.  +S2  k, 

1.5  k,  .k/ 

It  is  in  agreement  with  an  experimental  dependence  on  the  power  supply  source  energy  of  a 
total  substance  mass  evaporating  from  the  propellant  surface.  Indeed,  if  one  assumed  that  the 
propellant  surface  area  -  from  where  evaporation  takes  place  -remained  to  be  constant,  then. 


with  increase  in  the  energy  flux,  the  ratio  m/Wo  will  be  reduced  from  1.5  to  0.5  pg/J.  The 
energy  flux  density,  in  this  case,  is  increased  from  5.10"*  to  5.10®  J/cm^  s. 
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APPT,  the  considerable  energy  fraction  absorbed  by  the  propellant  surface  may  be  transferred 
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energy  flux  absorbed  by  propellant  surface  in  the  APPT  was  measured  with  a  small  size 
temperature  probes.  APPT  model  works  at  60  -300  J  bank  energy  with  average  mass  flow 
rate  of  order  10  ^^cm'^s'*.  Heat  flux  and  mass  flow  rate  measurements  were  analysed  on  the 
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power  consumption  (up  to  20-200  W).  Because  of  relative  simplicity  and  cheapness  of  SSC, 
simplicity  and  cheapness  of  PS  should  play  not  the  last  part  in  this  choice.  From  this  point  of 
view,  ablative  pulsed  plasma  thruster  (APPT)  having  many  engineering  and  operational 
advantages  is  the  most  perspective  option  for  SSC.  However,  low  thrust  efficiency  of  the 
widely  known  “traditional”  APPT  limits  its  real  capabilities  in  the  control.  Design, 
technologic  and  operational  simplicity  of  the  “traditional”  thruster  were  kept  in  the  APPT 
developed  during  the  few  previous  years.  Along  with  this,  new  thruster  is  characterized  by  a 
highly  efficient  operating  process  offering  thrust  efficiency  of  up  to  20-40%  within  the  stored 
energy  range  of  50-200  J.  This  paper  presents  the  development  and  refinement  results  for  the 
APPT  with  the  stored  energy  of  150  J  (APPT-150)  that  may  be  used  rather  efficiently  for  a 
control  within  the  power  consumption  range  of  20-200  W. 
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There  was  made  characterization  of  the  SPT-25  model  with  discharge  chamber  walls  made  of 
different  materials  such  as  the  Russian  BN-SiOi  (BGP)  type  ceramics,  AIN-BN  and  BN  types 
ceramics  produced  in  France  and  in  USA,  respectively.  The  SPT-25  model  has  external 
accelerating  channel  diameter  25mm  and  acceptable  performance  level  under  discharge 
powers  (1 00-200) W.  To  characterize  its  operation  with  discharge  chamber  walls  made  of 
different  ceramics  there  were  determined  the  voltage-current  characteristics  of  the  mentioned 
model  under  different  magnetic  fields  inside  the  accelerating  channel  and  different  mass  flow 
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rates  through  the  accelerating  channel.  There  were  determined  also  the  thrust  values  and  other 
output  parameters  such  as  the  thrust  efficiency  and  specific  impulse,  the  accelerated  ion  flow 
divergence  by  measurement  of  the  accelerated  ion  current  distributions  in  off-axis  angle  by 
RPA  probe,  some  plume  plasma  parameters  measured  by  cylindrical  electrostatic  probe  and 
discharge  voltage  and  current  oscillation  characteristics. 


264 


